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Summary 

The main objective of this study is to assess the specific capabilities of the defect energy parameter 
technique for global damage detection developed by Saleeb and coworkers. The feature extraction is the 
most important capability in any damage detection technique. Features are any parameters extracted from 
the processed measurement data in order to enhance damage detection. The damage feature extraction 
capability was studied extensively by analyzing various simulation results. The practical significance in 
structural health monitoring is that the detection at early stages of small-size defects is always desirable. 
The amount of changes in the structure’s response due to these small defects was determined to show the 
needed level of accuracy in the experimental methods. 

The arrangement of fine/extensive sensor network to measure required data for the detection is an 
“unlimited” ability, but there is a difficulty to place extensive number of sensors on a structure. Therefore, 
an investigation was conducted using the measurements of coarse sensor network. The white and the pink 
noises, which cover most of the frequency ranges that are typically encountered in the many measuring 
devices used (e.g., accelerometers, strain gauges, etc.) are added to the displacements to investigate the 
effect of noisy measurements in the detection technique. The noisy displacements and the noisy damage 
parameter values are used to study the signal feature reconstruction using wavelets. The enhancement of 
the feature extraction capability was successfully achieved by the wavelet theory. 


iii 
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Chapter I 
Introduction 

1.1 General 

Engineering structures are considered valuable assets to a nation due to their cost, complexity, and 
importance. These structures may deteriorate or damage during their service lives due to various reasons, 
that is, environment, loads, etc. The interest in the ability to monitor a structure and detect the damage at 
the earliest possible stage is common throughout the civil, mechanical, and aerospace engineering fields. 
Undetected damage may cause failure leading to loss of life and property. If the “state of health” of an 
important structure is constantly monitored, its maintenance cost can be kept to a reasonable minimum. 
The cost of repair is obviously less than that required to reconstruct the whole structure. As a result, 
structural health monitoring has received considerable attention in the literature. 

The increase in research activity regarding damage detection is the result of many factors. These 
factors can be generally categorized as spectacular failures resulting in loss of life, economic concerns, 
and recent technical advancements. Advancements in sensors and in the finite-element method have 
contributed to recent improvements in damage detection. Additional factors that have contributed to these 
improvements are the adaptation and advancements in experimental techniques such as modal testing and 
development of linear and nonlinear system identification methods. 

Damage is defined as changes introduced into a structure, either intentional or unintentional, which 
adversely affect the current or future performance of that structure. Damage may also be defined as any 
deviation in the structures original geometry or material properties that may cause undesirable stresses, 
displacements, or vibrations on the structure. Damage, such as a crack, in a structure reduces its stiffness. 
Similarly, the damping ratio will increase when damage progresses in the structure. The natural 
frequencies and the mode shapes are directly related to the stiffness of the structure. The loss of stiffness 
causes a drop in natural frequencies and a change in mode shapes. During the last decade, the use of 
dynamic system parameters such as natural frequencies, damping ratios, and mode shapes to detect 
damage has been studied intensively and many damage-detection techniques have been developed. 

Even though there has been significant research in the field of nondestructive damage identification, 
several factors remain to be tested to judge the potential of any damage-detection technique before it 
becomes a routine activity. Factors such as the limitation of measuring only a few modes in the 
experiments particularly in large structures, the limitation in the arrangement of an extensive sensor 
network, the experimental errors in the measurements, and the feature extraction capability to determine 
the location and the severity of the damage are to be investigated. Along with the above investigations, 
feature reconstruction methods used to extract the required signal from actual noisy data measurements 
need to be studied so as to enhance the feature extraction capabilities of the global damage detection 
techniques. 


1.2 Problem Statement 

Extensive research work over the past years using concepts of structural health monitoring have been 
directed towards the development of many damage-detection techniques. For the success of any 
technique, the following items should be studied: 

(1) The ability of the damage-detection scheme to withstand noisy measurements 

(2) Practicality of the density of the sensor network needed 

(3) Feature extraction and representation 

(4) Feature reconstruction with significant amounts of data/measurement noise 
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In the present study, the focus will be on the detail assessment of the four items listed above in 
connection with the specific defect energy parameter technique developed recently for global damage 
detection by Saleeb, et al. (refs. 67, 68, 80, and 85). In particular, extensive utilization of wavelet theories 
and graphical toolkits are made for the studies conducted under items (3) and (4) above. 


1.3 Objective of the Study 

The objective of this study is the assessment of the capabilities of the global damage detection 
technique developed by Saleeb and coworkers. This study presents the work of damage identification on 
plates when subjected to out-of-plane vibrations for a variety of cases. The assessment includes the 
validation of the technique by considering different types of failure/damage modes; the investigation on 
the density of the sensor network, where measurements on coarse sensor networks have particular 
importance; determination of the noise level that can be withstood by the detection technique; conduction 
of the false alarm tests; and the investigation of the damage feature extraction and representation for the 
above cases. An additional goal of this study is to reconstruct the damage feature from the noisy 
measurements by the utilization of the denoising tools, that is, wavelet transforms. 


1.4 Outline 

The introductory chapter is followed by background and review on literature in the areas of damage 
detection, data noise simulations, and the feature reconstruction using wavelet analysis in chapter II. 
This will follow the theoretical development on the proposed damage technique and wavelet theory in 
chapter III. Chapter IV describes the computational tools used to complete this work. The results and 
simulation case studies are presented in chapter V. Finally, a summary and conclusions are presented in 
chapter VI. 
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Chapter II 

Background and Literature Review 

2.1 Background 

Damage detection is a challenging problem that is under vigorous investigation by numerous research 
groups using a variety of analytical and experimental techniques. A significant amount of research work 
over the years has dealt with general topic of structural health monitoring. Health monitoring techniques 
may be classified as global or local. Global methods attempt to simultaneously assess the condition of the 
whole structure, whereas local methods provide information about a relatively small region of the system 
by using local measurements. Clearly, these approaches are complementary to each other, with the choice 
of methods being dependent on the scope of the problem at hand and nature of the sensor network on the 
structure. The damage detection procedure depends on the level of damage and deterioration of concern, 
the nature of the instruments, the spatial resolution of the sensors, the degree of the measurement noise, 
the configuration of topology of the test structures, complexity of the detection scheme, and the depth of 
knowledge concerning the failure modes of the structure. 

Another confounding factor is the fact that damage typically is a local phenomenon. Local response is 
captured by higher frequency modes whereas lower frequency modes tend to capture the global response 
of the structure and are less sensitive to local changes in a structure. From a testing standpoint it is more 
difficult to excite the higher frequency response of the structure, as more energy is required to produce 
measurable response at these higher frequencies than at the lower frequencies. These factors coupled with 
the loss of information resulting from the necessary reduction of time-history measurements to modal 
properties add difficulties to the process of vibration-based damage identification. 

The Non Destructive Evaluation (NDE) methods for detection of damage in structural systems have 
received increasing attention in the recent past for different classes of structural components using the 
signature analysis of the system response. While there are many approaches that have been investigated, 
or still being developed for signature based NDE of structure, the class of health-monitoring approaches 
that do not require detailed knowledge of the vulnerable parts of the structure, or of the failure mode of 
the structure, have a significant advantage in that they have the potential to cope with unforeseen failure 
patterns. Classes of NDE schemes that are less sensitive to the effects of initial assumptions are certainly 
of significant advantage. Considering the details of the NDE technique, focus is placed on the response 
signature-analysis, providing a global approach for fault diagnosis, as opposed to other local methods that 
require one to look directly at the location of suspected defects. 

The effect of damage on a structure can be classified as linear or nonlinear. A linear damage situation 
is defined as the case when the initially linear-elastic structure remains linear-elastic after damage. The 
changes in modal properties are a result of changes in the geometry and/or the material properties of the 
structure, but the structural response can still be modeled using linear equations of motion. Linear 
methods can be further classified as model-based and non-model-based. Model-based methods assume 
that the monitored structure responds in some predetermined manner that can be accurately discretized by 
finite element analysis, such as the response described by Euler-Bemoulli theory. 

Nonlinear damage is defined as the case when the initially linear-elastic structure behaves in a 
nonlinear manner after the damage has been introduced. One example of nonlinear damage is the 
formation of a fatigue crack that subsequently opens and closes under the normal operating vibration 
environment. Other examples include loose connections that rattle and nonlinear material behavior such 
as that exhibited by polymers. The majority of the studies reported in the technical literature address only 
the problem of linear damage detection. 

Another classification system for damage-identification methods defines four levels of damage 
identification, as follows: 

• Level 1 : Determination of the presence of damage in the structure. 

• Level 2: Level 1 plus determination of the geometric location of the damage. 
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• Level 3: Level 2 plus quantification of the severity of the damage. 

• Level 4: Level 3 plus prediction of the remaining service life of the structure. 

Level 1 is the most critical, whereas level 2 and level 3 would generally require additional more 
detailed measurements. Level 4 prediction is generally associated with fields of fracture mechanics, 
fatigue-life analysis, or structural design assessment. There are two alternative approaches for 
mathematical representation and implementation for these global detection techniques, i.e., system- 
identification and direct post processing of measurement data. The approaches mainly differ in the 
amount and type of information sets utilized. In particular, the system-identification approach is typically 
based on a complete analytical model that is fitted to the measured response 

Within the framework of system identification techniques, there is a distinction between methods that 
are used for continuous monitoring of structure performance and methods that are applicable to the 
detection of damage caused by extreme events. As an example, a system that uses continuous or 
intermittent accelerometer measurements from sensors mounted permanently to a bridge is different in 
terms of instrumentation and data acquisition requirements from a system that does not acquire data 
except during and immediately following an earthquake or a hurricane. It should be noted that primary 
distinction between these situations has to do with the sensors and data acquisition system requirements. 
Typically, the same types of analytical techniques can be applied to the data to determine the integrity of 
the structure. 

In addition to the intense computational demands, another main disadvantage of the system 
identification method is the need to treat the ‘inherent nonuniqueness’ caused by incompleteness of 
measured data with noise. Most recently, artificial neural networking has been used as an attempt to 
remedy some of these problems, but they still remain computationally intensive. 

On the other hand the direct detection scheme requires a key ingredient for the selection of the 
appropriate ‘measure’ or damage parameter, which is sufficiently sensitive to any slight perturbation in 
system properties. The whole implementation simply reduces to processing of raw experimental data, the 
objective being pattern-recognition, using the notion of comprised processed signatures, i.e., present 
(damage) vs. the base/reference (intact) states. Sharper resolution will be obtained if there are more 
distinct differences between two signatures. The presence of noisy data and consistency of distinctive 
patterns under different excitations (i.e., different vibration modes or different static load intensities) 
reflects on the robustness of NDE technique. It is this type of direct global detection methodology that is 
used in the later parts for the parametric and the assessment studies performed in this report. 

NDE technique should be robust enough to withstand certain level of unavoidable noise in practical 
measurements. If the technique is excessively sensitive to the noise levels then the effect of “true” 
damage will be completely “overshadowed” by the noise. With an eye towards reducing noise intensity, 
there has been development of different types of filters. The important factor in these so-called 
“denoising” schemes in the damage detection is retaining the significant, damage-produced features, of 
the signals. In the last two decades there has been extensive research going on the use of wavelet theory 
for this purpose. One of the applications of the wavelet theory is denoising. Wavelet theory is a useful 
tool to reduce the intensity of the noise in damage identification technique. Indeed, it is one of the main 
objectives in the present study to demonstrate the practical usefulness of applying the wavelet theories as 
a “denoising” tool in conjunction with a robust (global/direct) detection technique. 

With the above background in mind, we proceed next to summarize the pertinent literature. For 
convenience in the presentations, these are grouped into three separate parts; i.e., (a) Global Detection, (b) 
Data Noise Simulations/Filtering schemes as well as (c) applications of the wavelet theories in 
Denoising/Signal - Feature reconstructions. 
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2.2 Literature Review — An Overview 


The term damage refers to degradation or a failure of material. It can originate from diverse 
phenomena such as oxidation, carbonation, mechanical work, or any type or disintegration or weakening 
from aging or mechanical process. Within the framework of damage mechanics, only that which causes 
the loss of area, associated with change in local material properties such as Young’s modulus E, moment 
of inertia I, stiffness and flexibility and energy dissipation, is considered. When such changes occur there 
is a change of the entire physical system. This change leads to a change in vibration characteristics in 
physical space and also in material space. 

In physical space, when abrupt reductions in the cross section of a beam are considered, the properties 
have been changed, especially for the relationships between geometry and the centerline deflection 
curves. For free vibration behavior, such a beam results in a noticeable error in the natural frequencies 
because of the overestimation of the bending stiffness. It was found that resonant frequency and vibration 
amplitudes were considerably affected by the presence of cracks. Changes in the energy with respect to 
the damage are also important. Damage to engineering materials essentially results in a decrease of the 
free energy stored in a body with consequent degradation of the material stiffness. 

In material space, damage creates heterogeneity in a homogenous body, with consequent reduction of 
material properties including Young’s modulus, stiffness and moment of inertia. 


2.3 Literature Review — Detection Methods of System Identification Type 

Difficult challenges in formulating damage parameters possessing all these and other desirable 
attributes (e.g., applicability to different materials, multiple damage sites, various support condition, 
different vibration modes, etc.) are indicated by the numerous proposals made over the years (e.g., natural 
frequencies, mode shapes, influence flexibility coefficients, strain mode shape, curved mode shapes, as 
well as ratios, differences, fractions obtained form them). Basically, the performance of these measures 
was found to be heavily problem-dependent, with several conflicting conclusions often reached, when 
using the same measure under different conditions. 

Damage identification methods are based on the modification of structural model matrices such as 
mass, stiffness, and damping to reproduce as closely as possible the measured static or dynamic response 
data. These methods solve for the updated matrices (or perturbations to the nominal model that produce 
the updated matrices) by forming a constrained optimization problem based on the structural equations of 
motion, the nominal model, and the measured data. Comparisons of the updated matrices to the original 
correlated matrices provide an indication of damage and can be used to quantify the location and the 
extent of the damage. The methods use a common basic set of equations, and the differences in the 
various algorithms can be classified as follows: 

1. Objective function to be minimized 

2. Constraints placed on the problem 

3. Numerical scheme used to implement the optimization. 

Chen and Garba (ref. 1) presented a method for minimizing the norm of the model property 
perturbations with a zero modal force error constraint. They also enforce a connectivity constraint to 
impose a known set of load paths onto the allowable perturbations. The updates are thus obtained at the 
element parameter level, rather than at the matrix level. This method is demonstrated on a truss FEM. 

McGowan, et al.(ref. 2) report ongoing research that examines stiffness matrix adjustment algorithms 
for application to damage identification. Based on the measured mode shape information from sensor 
locations that are typically fewer than the DOF in an analytical model, mode shape expansion algorithms 
are employed to extrapolate the measured mode shapes such that they can be compared with analytical 
model results. These results are used to update the stiffness matrix while maintaining the connectivity and 
the sparsity of the original matrix. 
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Smith (ref. 3) presented an iterative approach to the optimal update problem that enforces the sparsity 
of the matrix at each iteration cycle. Multiplying each entry in the stiffness update by either one or zero 
enforces the sparsity pattern. Kim and Bartkowicz (ref. 4) investigated damage detection capabilities with 
respect to various matrix update methods, model reduction methods, mode shape expansion methods, 
number of damaged elements, number of sensors, number of modes, and levels of noise. The authors 
developed a hybrid model reduction/eigenvector expansion approach to match the order of the 
undamaged analytical model and the damage test mode shapes in the matrix update. They also introduced 
a more realistic noise level into frequencies and mode shapes for numerical simulation. From both 
numerical and experimental studies, the authors showed that the number of sensors is the most critical 
parameter for damage detection, followed by the number of measured modes. Lindner, et al. (ref. 5) 
presented an optical update technique that formulates an over determined system for a set of damage 
parameters representing reductions in the extensional stiffness values for each member. The value 
represents the amount of stiffness reduction in that member. 

Lie (ref. 6) presented an optimal update technique for computing the elemental stiffness and mass 
parameters for a truss structure from measured modal frequencies and mode shapes. The method 
minimizes the norm of the modal force error. The author demonstrates that if sufficient modal data is 
variable, the elemental properties can be directly computed using the measured modal frequencies, 
measured mode shapes, and two matrices which represent the elemental orientations in space and the 
global connectivity of the truss. In this case, the solution for the elemental properties is shown to be 
unique and globally minimal. The method is used to locate a damaged member in a FEM of a truss using 
the first four measured modes in sets of three at a time. 

Zimmrman and Kaouk (ref. 7) presented the basic minimum rank perturbation theory (MRPT) 
algorithm. A nonzero entry in the damage vector is interpreted as an indication of the location of the 
damage. The resulting perturbation has the same rank as the number of modes used to compute the modal 
force error. It is demonstrated that the MRPT algorithm preserves the rigid body modes of the structure 
and the effects of measurement and expansion errors in the mode shapes are demonstrated and discussed. 

Kaouk and Zimmerman (ref. 8) further developed this algorithm and demonstrated how perturbations 
to two of the property matrices can be estimated simultaneously by using complex conjugates of the 
model force error equation. The method is demonstrated numerically for a truss with assumed 
proportional damping. Also, the technique is used experimentally to locate a lumped mass attached to a 
cantilevered beam. 

Kaouk and Zimmerman (ref. 9) extended the MRPT algorithm to estimate mass, stiffness, and 
proportional damping perturbation matrices simultaneously. The computation of these individual 
perturbation matrices is accomplished by exploiting the cross-orthogonally conditions of the measured 
mode shapes with respect to the damage property matrices. The authors examined the results by 
computing a cumulative damage vector. 

Kaouk and Zimmerman (ref. 10) presented a technique that can be used to implement the MRPT 
algorithm with no original FEM. The technique involves using a baseline data set to correlate an assumed 
mass and stiffness matrix, so that the resulting updates can be used as the undamaged property matrices. 

Zimmerman and Simmermacher (refs. 1 1 and 12) computed the stiffness perturbation resulting from 
multiple static load and vibration tests. This technique is proposed partially as a method for circumventing 
the mismatch in the number of models between test and FEM. They applied this technique to a FEM of a 
structure similar to a NASA test article. They also presented two techniques for overcoming the rank 
deficiency that exists in the residual vectors when the results of one static or modal test are linear 
combinations of the results of previous tests. 

Ricles (ref. 13) presented a methodology for sensitivity-based matrix update, which take into account 
variations in system mass and stiffness, center of mass locations, changes in natural frequency and mode 
shapes, and statistical confidence factors for structural parameters and experimental instrumentation. The 
method uses a hybrid analytical/experimental sensitivity matrix, where the modal parameter sensitivities 
are computed from the experimental data, and the matrix sensitivities are computed from the analytical 
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model. This method is further developed and applied to more numerical examples by Ricles and 
Kosmatka. 

Sanayei and Onipede (ref. 14) presented a technique for updating the stiffness parameters of a FEM 
using results of a static load-displacements test. A sensitivity -based, element-level parameter update 
scheme is used to minimize the error between the applied forces and forces produced by applying the 
measured displacements to the model stiffness matrix. The sensitivity matrix is computed analytically. 

The structural DOF are partitioned such that the locations of the applied loads and the locations of the 
measured displacements are completely independent. The technique is demonstrated on two FEM 
examples. 

In a related paper, Sanayei, et al. (ref. 15) examined the sensitivity of the previous algorithm to noisy 
measurements. The influence of the selected measurement DOF set on the error in the identified 
parameters is studied. A heuristic method is proposed that recursively eliminates the measurement DOF 
that the elemental stiffness parameters are the most sensitive. In this manner, the full FEM DOF set is 
reduced to a manageable size while preserving the ability to identify the structural stiffness parameters. 

Hemez and Farhat (ref. 16) presented a sensitivity -based matrix update procedure that formulates the 
sensitivities at the element level. This has the advantage of being computationally more efficient than 
forming the sensitivities at the global matrix level. It also allows the analysis to “focus” on damage in 
specific members. 

Zimmer and Kaouk (ref. 17) implemented an eigenstructure assignment technique for damage 
detection. They included algorithms to improve the assignability of the mode shapes and preserve sparsity 
in the updated model. They applied their technique to the identification of the elastic modulus of a 
cantilevered beam. 

Linder and Goff (ref. 1 8) defined damage coefficients for each structure member. They then used an 
eigenstructure assignment technique to solve for the damage coefficient for each member. They applied 
this technique to detect simulated damage in a 10-bay truss FEM. 

Schulz, et al. (ref. 19) presented a technique similar to eigenstructure assignment known as “FRF 
assignment”. The authors formulated the problem as a liner solution for element-level stiffness and mass 
perturbation factors. They pointed out that using FRF measurements directly to solve the problem is more 
straightforward than extracting mode shapes. They used measured mobility functions (FRFs from velocity 
measurements) to obtain higher numerical accuracy, since the velocity response is flatter over the entire 
spectrum than either the displacement or acceleration response. The technique is applied to an FEM of a 
bridge structure. 

Baruh and Ratan (ref. 20) used the residual modal force as an indicator of damage location. They 
separated the residual modal force into the effects of identification error in the measurements, modeling 
error in the original structural model, and modal force error resulting from structural damage. They 
examined the sensitivity of the damage location solution to errors in the original structural model and to 
inaccuracies in the modal identification procedure. 

Kim and Bartowicz (ref. 21) and Kim, et al.(ref. 22) presented a two-step damage-detection procedure 
for large structures with limited instrumentation. The first step uses optimal matrix update to identify the 
region of the structure where damage has occurred. The second step is a sensitivity-based method, which 
locates the specific structural element where damage has occurred. The first advantage of this approach 
lies in the computational efficiency of the optimal update method in locating which structural parameters 
are potentially erroneous. The second advantage lies in the small number of parameters updatation by the 
sensitivity-based technique. 

Li and Smith (ref. 23) presented a hybrid model update technique for damage identification that uses 
a combination of the sensitivity and optimal-update approaches. This method constrains the stiffness 
matrix perturbation to preserve the connectivity of the FEM, and the solution minimizes the magnitude of 
the vector of perturbations to the elemental stiffness parameters. 

Hung-Liang, et al. (ref. 24) presented a nondestructive evaluation method to identify the structural 
stiffness of ceramic candle filters. All filters were subjected to an excitation force, and the response was 
picked up by an accelerometer in a free-free boundary condition. The frequency response function and 
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vibration mode shapes of each filter were evaluated. The results are also estimations of the overall 
bending stiffness values for four different types of candle filters. The used filters showed stiffness 
degradation. The location and the size of the damaged section were identified using the measured modal 
strain energy. 

Dynamic bending stiffness was used by J. Maeck et al. (ref. 25) to detect damage. They discussed 
different techniques and derived from experimentally determined modal characteristics of a reinforced 
concrete beam from its dynamic bending stiffness. The degradation of stiffness, due to cracking of the 
reinforced concrete, gives information on the position and severity of the damage that has occurred. 

Yuen (ref. 26) presented a systematic study of the relationship between damage location and size, and 
the change in the eigenvalues and eigenvectors of a cantilever beam. Damage was modeled as a modulus 
reduction in an element of the beam. The eigen parameters were studied, i.e., translation and rotation. 

Both showed a sudden change at the damage region along the beam coordinate. However, the rotation 
eigen parameter was not detected for higher modes. 

Ren and De Roeck (refs. 27 and 28) proposed a damage identification technique based on a change in 
frequencies and mode shapes of vibration, for predicting damage location and severity. The method is 
applied at an element level with a conventional finite element model. The element damage equations have 
been established through the eigen- value equations that characterize the dynamic behavior. The influence 
of noise was also shown and they verified their method by a number of damage scenarios for simulated 
beams and found the exact location and the severity of damage. They demonstrated that multiplying the 
damage eigen-value equations with the undamaged or damaged mode shapes provides more equations 
and guarantees the damage localization. 

Stubbs and Topole (ref. 29) proposed a formulation that localizes and determines the size changes in 
the stiffness of the structure. Generally such changes are a reduction in stiffness and are associated with 
some type of structural damage. Serious damage will change the stiffness locally and globally. Thus a 
reduction in stiffness is generally interpreted as damage. However, reductions in stiffness do not 
necessarily relate to damage. Therefore, the algorithm is a conservative method to determine potential 
locations of damage. 

Gawronski and Sawicki (ref. 30) used modal and sensor norms to determine damage locations in 
flexible structures. It provided information about the impact of the damage on the natural modes of the 
damage structures. As the norm is determined from the system natural frequencies, modal damping ratios, 
and the input and output gains; they depend neither on the input time history nor the actual system 
deformation. 

From fracture mechanics concepts, Rizo et al.(ref. 31) developed a spectral method to identify cracks 
in various structures. Crack depth was related to the change of natural frequencies, but this method lacks 
accuracy for small cracks, and it is not quite clear that this approach can be generalized to complex 
structures. 

F. Vestroni and D. Capecchi (ref. 32) found damage by frequency measurement. A linear behavior 
was assumed, before and after the damage. The method was described and used when frequencies are the 
observed quantities. The procedure is generalized by assuming finite-element interpretative models and 
an automatic algorithm of modal updating, which is used to determine the best stiffness distribution for an 
assigned location of damage. A minimum amount of frequencies is necessary to obtain a unique solution. 
This is important, because the problem often over determined. Quantity of measured data is important to 
reach an acceptable solution. 


2.4 Literature Review — Detection Methods of Direct/Pattern Recognition Global Type 

One approach to structural health monitoring is to look for changes in a “signature” of the structure 
that is related to its dynamic characteristics. The observed changes in the structure, for example modal 
parameters, are compared to a database of possible changes and the most likely change is selected for 
detecting damage and locating its position. It is first noted that for beams, plates, and shells there is a 
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direct relationship between curvature and bending strain. Some researchers discuss the practical issues of 
measuring strain directly or computing it from displacements or accelerations. 

Stubbs, et al.(ref. 33) presented a method based on the decrease in modal strain energy between two 
structural DOF, as defined by the curvature of the measured mode shapes. Topole and Stubbs (refs. 34 
and 35) examined the feasibility of using a limited set of modal parameters for structural damage 
detection. They studied several dynamic parameters for damage detection using full scale modal testing. 

A probable failure due to a large fatigue crack was simulated by unfastening a set of high-strength bolts in 
a splice connection of a steel highway bridge. Experimental modal testing was performed for intact case 
as well as the cracked case. Results indicate the presence of detectable changes in some of the response 
data to a simulated physical failure. Non-parametric information, i.e., time records, frequency spectra, 
transfer functions as well as parametric information, i.e., modal frequencies and mode shapes had been 
examined. The Modal Assurance Criteria (MAC) and modal frequencies can detect the damage in higher 
modes; otherwise the modal frequencies were not sensitive. In more recent publications, Stubbs and Kim 
(ref. 36) examined the feasibility of localizing damage using this technique without baseline modal 
parameters. 

Pandey, et al.(ref. 37) developed curvature-mode shapes in which the absolute change is located in 
the region of damage, and hence can be used to detect damage in a structure. The changes in the curvature 
mode shape increases with increasing damage. The difference in modal curvature between the intact and 
the damaged beam showed not only a high peak at the fault position, but also some small peaks at 
different undamaged locations for higher modes. This can cause confusion to the analyst in a practical 
application in which one does not know in advance the location of faults. Also, it is not sensitive for small 
damages. 

Chance, et al. (ref. 38) found that numerically calculating curvature from mode shapes resulted in 
unacceptable errors. They used measured strains instead to measure curvature directly, which 
dramatically improved results. 

One feasible measurement property is change of energy with respect to material damage. Damage to 
engineering materials results in a decrease of the material stiffness. DisPasquqle, et al. (ref. 39) found that 
the parameter based global damage indices could be related to locate damage variables through operations 
that average over the body volume. 

Stubbs and Osegueda (refs. 40 and 41) evaluated non-destructive damage detection. The damage was 
modeled by reduced modulus in the element. A finite element model was used to establish the sensitivity 
matrix. Laboratory experiments were performed for supporting the concept. As far as the damage location 
is concerned, false predictions still occurred. Moreover, from the experimental results, small damage 
levels cannot be detected. 

West (ref. 42) presented what is possibly the first systematic use of mode shape information for the 
location of structural damage without the use of a prior FEM. The author uses the Modal Assurance 
Criteria (MAC) to determine the level of correlation between modes from the test of an undamaged Space 
Shuttle Orbiter body flap and the modes from the test of the flap after it has been exposed to acoustic 
loading. The mode shapes are partitioned using various schemes, and the change in MAC across the 
different partitioning techniques is used to localize the structural damage. 

Mayes (ref. 43) presented a method for model error localization based on mode shape changes known 
as structural translational and rotational error checking (STRECH). By taking ratio of relative modal 
displacements, STRECH assess the accuracy of the structural stiffness between two different structural 
degrees of freedom (DOF). STRECH can be applied to compare the results of test with an original FEM 
or to compare the results of two tests. 

Ratcliffe (ref. 44) presented a technique for locating damage in a beam that uses a finite difference 
approximation of Laplacian operator on mode shape data. Cobb and Liebst (ref. 45) presented a method 
for prioritizing sensor locations for structural damage identification based on an eigenvector sensitivity 
analysis. 

Modal curvature method was used by Wahab, et al.(ref. 46) to detect damage in a pre- stressed 
concrete bridge. To establishment their method they used simulated data from simply supported and 
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continuous beams containing damage parts at different locations. A damage indicator called a “curvature 
damage factor” was introduced which is the difference in curvature mode shapes for all modes and could 
be summarized by one number for each measured point. For several damage locations in the structure, all 
modes should be carefully examined. The lower modes are, in general, more accurate than the higher 
modes. When more than one fault exists in the structure, it is not possible to locate damage in all positions 
from the result of only one mode. 

Wavelet-based approach for Structural Health Monitoring (SHM) and damage was used by Z. Hou, et 
al. (ref. 47). The method is applied to simulated data generated from a simple structural model subjected 
to a harmonic excitation. Spikes in details of the wavelet decomposition may detect changes in system 
stiffness, and the locations of these spikes indicate the moment when structural damage occurred. 

J.C. Hong, et al. (ref. 48) applied wavelet transforms for the damage identification in a structural 
member. They showed the effectiveness of the wavelet transforms by means of its capability to estimate 
the Lipshitz exponent. A model beam which has a defect represented by an abrupt thickness change is 
taken for the study. They derived bound of the Lipschitz exponent using the simplest Euler beam theory. 
Based on the estimated Lipschitz exponent, they suggested an optimal wavelet for the estimation of the 
Lipschitz exponent. The Continuous Wavelet Transform (CWT) by a Mexican hat wavelet having two 
vanishing moments is utilized. CWT is applied on the mode shape of the damaged beam. The modulus 
maximum lying inside of the cone of influence is clearly seen in the contour plot of CWT. It is also able 
to detect the location of the damage. 

Y.J. Yan, et al.(ref. 49) presented structural damage feature index using dynamic response of the 
structure. They used honeycomb sandwich plate with free-free boundary conditions. Dynamic responses 
are measured using piezo-patch sensor. Wavelet transforms are utilized to calculate the damage index. 
Dynamic response signal of the plate is decomposed into 5 th layer of wavelet transform, and 16 sub- 
signals are obtained. The relative energy variation in each sub-signal is calculated. This approach is 
evaluated by comparing simulated and experimental data. The error between simulated and experimental 
data energy variation is less than 10 percent with a crack length more than 3 percent of the length of the 
plate. This method can able to tell whether the damage is exist or not, but it can not able to tell the 
location and severity of the damage. 


2.5 Literature Review - Data Noise Simulation and Filtering Schemes 

When the measured data is noise free, an element can be viewed as damaged if its estimated 
parameter is different from the baseline value. This simple assessment is complicated by the presence of 
measurement noise. The noise in the measurements would cause the estimated parameter for an element 
to be different from the baseline value even if there is no damage at all. 

Ikumasa Yoshida and Tadanobu Sato (ref. 50) suggested from the standpoint of damage detection, 
non-Gaussian non-white noises might be preferable, because the damage tends to be concentrated on a 
specific part of a structure. From this consideration, exclusive noise which is non-Gaussian and non-white 
is proposed. They explained from observations of actual damages, proper process noise for damage 
detection should have following nature, i) concentration in space, ii) correlation in time domain. Proposed 
exclusive noise has both nature i) and ii). The exclusive noise is given to only one element at each step, so 
that it is consistent with the damage concentration nature. At the next step, the element is given the noise 
is selected with probability Ps , otherwise an element is selected randomly. Consequently, the noise has 
correlation in time domain. 

Taeho Charles Jo (ref. 51) presented on MLP training with noise optimized by genetic algorithm. He 
explained the noise addition to the training patterns as follows. The input vector x t is represented x t = \x ilf 

x i2 , , x im \. The dimension of input vector is m. The noise injected into the input vector x t of training 

pattern and 8/ is represented as 8/ = [ 8^8^ , z im ]. The element of 8/ is the random number 
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based on Gaussian distribution. Therefore, the noise is injected into the input vector like x\ = 

[ x i\ + £i\, x ii + £/2, , x im + ]. In this case, the input vector with noise is x / . 

Xu and Liu (ref. 52) developed a combined optimization technique of using an improved p GA and 
local optimizer to solve the optimization problem so as to obtain the flaw parameters defining flaw 
configuration. To investigate the effect of noise, they introduced white noise r(x) intentionally in the 


n 

simulated measurements for the surface displacement responses. T(x) is described as ^ r(x z ) = 0, 

i = 1 


1 

that the improved p GA can bear relatively weak noise disturbances of r| = 5 percent. To effectively 
suppress the effect of noise, they suggested to filter the noisy displacement responses before they are used 
in the objective function. After filtering the noisy data they can able to go up to r| =15 percent. 

Cacciola, Impollonia, and Muscolino (ref. 53) adopted Monte Carlo Method for the crake detection in 
a damaged beam. They explained the analytical model of the cracked beam and considered the real-life 
test data which will contain the noise. In their paper the stochastic dynamic analysis of a cracked 
cantilever beam under white noise is addresses in the time domain. 

Prashant Pawar and Ranjan Ganguli (ref. 54) proposed a genetic fuzzy system for damage detection 
in beams and helicopter rotor blades. They added noise in data to simulate the uncertainty present in 
experimental measurements and the modeling process. Given a computed frequency measurement delta 
Aco , random number u in the interval [-1, 1] and a noise lever parameter a , the noisy simulated data is 
given as A(jd noisx = A<x>f em + u a . The parameter a defines maximum variance between the computed 


yrfe)ffe - t) = 2D5(x) and D = r \ 


A/2 


i = 1 


x i 


, where r| is the noise amplitude. It can be seen 


value of Aco and simulated measured value A(ti noisx which is a simulation of practical measurement. 
A<x>f em is measurement delta obtained by finite element analysis. They conclude that the genetic fuzzy 
system gives an average success rate in damage detection of 99.81 percent, even with the noise levels of 
0.20 in the data. For noise levels of 0.10 and 0.15, the average success rate is about 100 percent. The 
genetic fuzzy system is therefore very robust in the presence of significant noise inherence in mechanical 
system measurements. 

Among ecologists, there has been a growing recognition of the importance of long-term correlations 
in environmental time series. John Halley (ref. 55) described that the family of 1 If noise-fluctuations 
defined in terms of the different timescales present- is a useful approach to this problem. White and 
random walk, the two currently favored descriptions of environmental fluctuations, lie at extreme ends of 
this family of process. Ecologists expect both rare and common events to be important. 1 If noise is a way 
describing these kinds of events. It is associated with an altogether slower decline in correlation. The 


extended family of 1 If noise is characterized by power-law spectra of forms: Spectral density, S(f) 



where 0 < y < 2. S(f) is a decomposition of the noise signal into various component frequencies per unit 
frequency. By analogy with light, white noise is noise, whose spectral density is flat, containing equal 
amounts of all frequencies. The 1 If family takes its name from that member for which y « 1, which is 
often called ‘pink noise’. It shows no preference for short or long timescale disturbances. He suggested 
that pink 1 If noise, which lies midway between white noise and random walk, might be the best null 
model of environmental variation. 

1 If noise is widely present in electrical components. The presence of the excess low-frequency noise 
is particularly troublesome for microelectronics reliability, because its amplitude is much larger than that 
of other limiting noises. Keiji Takagi (ref. 56) provided experimental evidence of low-frequency noise in 
granular resistors. He measured the noise level with a spectrum or an FFT analyzer and a noise intensity 
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measuring system and noticed that all noise spectra are of the 1 If type in the low-frequency range. He 
explained a model of Hf noise spectrum generation in granular structures. So there can be 1 If noise in the 
experiments using electronic devices. 

Widespread occurrence of signals exhibiting power spectral density with 1 If noise behavior suggests 
that a general mathematical explanation of such an effect might exits. Mathematical algorithms and 
models for the generation of the processes with 1 If noise cannot, as a rule, be solved analytically and they 
do not reveal the origin as well as the necessary and sufficient conditions for the appearance of 1 If type 
fluctuations. Kaulakys and Meskauskas (ref. 57) proposed a simple analytical solvable model of 1 If noise. 
The models reveal main features, parameter dependencies and possible origin of 1 If noise, i.e., random 
increments of the time intervals between the pulses. The conclusion that l/ / 5 noise with 8 « 1 may result 
from clustering of the signal pulses, particles or elementary events can be drawn from the analysis of the 
simple, exactly solvable models. The mechanism of clustering depends on the system. 

Rangarajan (ref. 58) explained about random noise, structured noise, and physiological interference. 
Random noise: it refers to an interference that arises from a random process such as thermal noise in 
electrical devices. A random process is characterized by the probability density function (PDF) 
representing the probabilities of occurrence of all possible values of a random variable. It is common to 
assume the mean of random noise process to be zero. In most cases the noise is additive, y(t ) = x(t) + n(t). 
There are cases where the noise can be multiplicative, y(t ) = x(t) * n(t) where y(t) is measures signal, x(t) 
is original signal, and n(t) is the noise. Structured noise: Power-line interference at 50 or 60Hz is an 
example of structured noise: the typical waveform of the interference is known in advance. It should, 
however, be noted that the phase of the interfering waveform will not usually be known. Physiological 
interference: The human body is a complex conglomeration of several systems and processes. Several 
physiological processes could be active at a given instant of time, each one producing many signals of 
different types. Physiological interference may not be characterized by any special waveform or spectra 
content, and is typically dynamic and non-stationary. 

A fast growing variety of computer applications require use of random numbers. An immediate 
example can be found in a source of noise superimposed on the signals in studies of digital signal 
processing. Such random numbers are generated by iterative calculations based on a numerical formula. 
The sequences, generated based on the multiplicative congruential algorithm and its modifications that are 
often used in the functions generating random numbers for PCs are characterized periodically. These 
random numbers are not completely random. In view of the above mentioned facts, Shunsuke Kishimoto 
and Masuo Fukue (ref. 59) decided to develop so called “physical random digits,” that is random digits 
generated by hardware and intended for safe use in PCs. The idea proposed by them consists in placement 
of all the hardware necessary for generation of physical random digits on a single IC chip. For this 
purpose, they used a Zener diode (semiconductor diode) that can be incorporated in the IC chip as a 
source of noise, and attempted to generate sequences of random digits based on the randomness of the 
short noise generated by the diode within time intervals of a fixed duration. This new method was verified 
by generating 1.2xl0 9 decimal random integers and by testing such characteristics of the generated 
random number sequences as their uniformity and independence by statistical methods. 

In control experiment, degraded images are obtained by adding noise to the original image and then 
they are restored by the restoration technique. In most cases, it is assumed that the noise is uncorrelated to 
the original image. Shan Suthaharan and Ray (ref. 60) introduced a method to generate a noise (called the 
first noise) such that average cross power spectrum between the first noise and the original image is zero. 
This method is extended to generate the second noise such that the average cross power spectrum between 
the second noise and the original image and the average cross power spectrum between the first and 
second noise are zero. Let / represent the original image. Suppose we generated, arbitrarily, an initial 
noise wl. Because of the very nature of its being generated arbitrarily, w 1 could be correlated with f To 
make the noise uncorrelated with f they define first noise nl as n\ = wl - ( Cov (wl, f)/Var(f))*f 
Suppose a second initial noise w2 is generated arbitrarily. Because of the very nature of its being 
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generated arbitrarily, w2 could be correlated with / and nl . To make the noise uncorrelated they define 
second noise nl as nl = w2- ( Cov (w2, nl)/Var(nl))*nl - ( Cov (wl, f)/Var(f))*f 

A single set of noise-polluted data can lead to a biased result. The effect of the noise on the topology 
of the objective function is evident in the sensitivity of parameter estimates when noisy data are subjected 
to random perturbation. Pothisiri and Hjelmstad (ref. 61) adopted a data perturbation scheme to create a 
sample of artificial data sets for parameter estimation at each stage of the damage localization process. 

They explained the solution multiplicity for the noisy data. To wit, the j th component of perturbed 

~ a 

eigenvector from the j th component of the z th measured eigen vector as §y = 4> z y (1 + r| z y ) , where r| y is a 

uniform random variate in the range [- a , a\. The amplitude a must be specified to account for the 
sensitivity of parameter estimate due to noise. Once the multiple solutions have been located using the 
random starting points, the data perturbation scheme can be applied to obtain the sensitivity information 
for each of the solutions by using the individual solution as a fixed starting point for each perturbed data 
set. 


2.6 Literature Review — Application of Wavelet Theory 
for Denoising and Signal Feature Reconstruction 

In all the experiments, presence of noise is a major issue. Denoising of the signals and images has 
taken major role in the current research activities. In the last two decades, vast research is going on 
wavelet theory and their applications. One of the applications of the wavelets is denoising of the signals 
and images. The wavelet transform has proven to be very successful in making signal and noise 
components of the signal distinct. The signal is transformed into some domain where the noise 
component is more easily identified, thresholding operation is then applied to remove the noise, and 
finally the transformation is inverted to reconstruct a (hopefully) noise-free signal. 

In the recent years there has been a fair amount of research on wavelet thresholding and threshold 
selection for signal de-noising, because wavelet provides an appropriate basis for separating noisy signal 
from the signal. Lakhwinder Kaur, Savita Gupta, and Chauhan (ref. 62) proposed a near optimal threshold 
estimation technique for image denoising which is subband dependent i.e., the parameter for computing 
the threshold are estimated from the observed data, one set for each subband. The proposed method is 
called “NormalShrink”. The threshold value (TV) is computed by T N = |3 a 2 / a y . P is scale parameter, 

a 2 is variance and a y is standard deviation of the subband. The image denoising algorithm is described 

and experiments are conducted to assess the performance of NormalShrink. The results show that it 
removes noise significantly. 

Under the framework of VC-theory, wavelet thresholding amounts to ordering of wavelet coefficients 
according to their relevance to accurate function estimation, followed by discarding insignificant 
coefficients. Existing wavelet thresholding methods specify an ordering based on the coefficient 
magnitude, and use threshold(s) derived under Gaussian noise assumption and asymptotic settings. In 
contrast, Shi Zhong and Vladimir Cherkassky (ref. 63) proposed a wavelet threshold approach, uses 
orderings based on the coefficient magnitude, which better reflects statistical properties of natural images, 
and VC-based thresholding developed for finite sample settings under very general noise assumptions. 
They interpreted image denoising as a special case of signal estimation problem and propose a model 
selection based denoising method under the framework of VC theory, which was developed for 
estimating data dependencies from the finite samples. VC-theory provides a framework for model 
selection called structural risk minimization (SRM). Under SRM, a set of possible models ordered 
according to their complexity. Level dependent thresholding has been proposed to improve the 
performance of wavelet thresholding method. Instead of using a global threshold, level-dependent 
thresholding uses a group of thresholds, one for each scale level. 
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Microarray imaging is a recent cutting-edge technology in bioinformatics which can monitor 
thousands of genes simultaneously. This technology can monitor thousands of DNA sequences in a high 
density array on a glass. The noise source in microarray imaging originates from different sources during 
the course of experiment, such as photon noise, electron noise, laser light reflection, dust on the slide, and 
so on. Hence, it is crucial to denoise the resultant image within this process. Wang, Robert Istepanian, and 
Yong Hua Song (ref. 64) proposed a new approach on wavelet theory to provide an enhanced approach 
for eliminating such noise source and ensure better gene expression. They apply the SWT (Stationary 
Wavelet Theory) method to preprocess the microarray images for removing the random noise. In this 
method at each level, when the high-pass and low-pass filters are applied to the data, the two new 
sequences have the same length as the original sequences. To do this, the original data is not decimated. 
However, the filters at each level are modified by padding them out with zeros. They apply the 
biorthogonal wavelet procedure to decompose the image using Matlab. SureShrink thresholding 
algorithm is applied to the subimages. The results show that SWT provides a better performance over 
traditional wavelet transform method. 

Ultrasonic is the perfect means for morphological investigation of the neonatal brain. A problem, 
common to all medical ultrasound images, is the presence of speckle noise, which not only complicates 
the visual interpretation of images, but also the quantitative measurements. Suppression of speckle in 
noise is necessary to get reliable measurements. Ivana Duskunovic, et al.(ref. 65) compared two noise 
removal techniques, based on wavelet decomposition, applied to speckle images. The first method applies 
Bayesian shrinkage of the wavelet coefficients using MRF model to express the prior knowledge about 
the spatial clustering of the coefficients. This method works as follows: Take the discrete wavelet 
transform of the observations, apply a simple nonlinearity (shrink towards zero) to each wavelet 
coefficient, and compute signal estimates by applying the inverse transform to the transformed 
coefficients. The second method proposed is based on the spatial analysis of detail images and 
determining the position of real edges and removing false edges. Values of the pixels, for which they 
determine that do not belong to an edge, are set to zero. Specifically, their method operates on detail 
images. Each technique has its advantages. The first technique provides significant noise suppression, 
which results in very smooth areas uniformly corrupted by noise. The second technique does not suppress 
as well as the first but preserves sharpness better. The techniques proposed yield comparable results. 

Gaussian smoothing is a simple denoising method and it is widely used in computer vision 
algorithms. However, Gaussian smoothing destroys the structure of image data. The orthonormal basis 
property of wavelet brings the possibility of denoising without destroying the structure of input data. Qi 
Li, et al. (ref. 66) investigated how useful the wavelet denoising technique is in improving the accuracy of 
computing fundamental matrices. They applied a methodology called waveShrink, which has wide 
applications and implemented in commercial software, e.g., wavelet toolbox of Matlab. Determining 
threshold is key issue in waveShrink denoising. They applied minimax threshold technique to denoise the 
input data. The experimental results show that wavelet denoising can help improving accuracy of 
fundamental matrix on a large scale of images while it is useless for stereo images containing randomized 
information or highly regular texture information. 


2.7 Concluding remarks 

A great deal of research in the past thirty years has been aimed at establishing an effective method for 
health monitoring in civil, mechanical, and aerospace structures. The objective is to determine the 
existence, location, and degree of damage in a structure. The development of a successful technology for 
structural health monitoring has enormous potential for application in evaluation of offshore structures 
and bridges subject to fatigue, corrosion, impact and earthquakes as well as buildings and aerospace 
structures subject to severe loads or structural deterioration. A variety of methods for evaluating damage 
in structural systems have emerged and evolved. All of these methods require a parameter estimation 
algorithm to drive them; i.e., the selection of an appropriate “measure” or suitable perturbation in system 
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properties. To this end, recent work by Saleeb, et al. (ref. 67) has been directed towards the development 
of appropriate global indices of this type, based on a more fundamental approach in structural mechanics. 

The novelty of this approach stems from the use of an alternative formulation in material (vs. physics) 
space. These dual balance laws are revealing in that the resulting force/source terms are directly and 
explicitly driven by increased heterogeneity due to deterioration. It thus provides an ideal candidate for 
the damage detection parameter. Such a damage parameter is easily computed from the measured raw 
data (e.g., strains, deflections, and rotations). A large number of numerical simulations and comparison to 
data have clearly demonstrated the power of this formulation under different test conditions (refs. 68 
and 69). 

Most of the experiments contain some amount of noise in the measured data. The damage 
identification formulation should be robust enough to sustain the noise to identify the damage. For the 
past few decades there has been decent amount of research in the field of development of the filters to 
reduce the noise in signal processing and image processing. Few filtering techniques are developed for 
special applications. The important feature of the filter is to retain the significant details. For the last two 
decade there is great deal of research is going on wavelet theory. One of the applications of the wavelets 
is denoising the data. It has the feature to retain the significant details. Plenty of wavelets are developed 
by the researchers, but all wavelets will not fit to particular data. Hence selecting the suitable wavelet for 
the study takes major role. The focus of this study is to provide further validation for the feasibility and 
effectiveness of this new algorithm under realistic conditions. Algorithm is evaluated with noisy data. 
Wavelet theory is adopted to reduce the intensity of the noise in the simulated data. 
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Chapter III 

Theoretical Developments 


The theoretical background necessary for this study is described in this chapter. This chapter is 
divided into two sections; Damage Identification Technique and Wavelet Theory. The damage 
identification technique that will be presented is utilized for the study of damage detection as applied to 
plate like structures and a unique application of wavelet theory is used to reduce the noise intensity of the 
measured data required in the detection scheme. 


3.1 Damage Identification Technique 

3.1.1 Introduction. — This section deals with the description of the underlying mathematical structure 
of the damage index used in the detection scheme and various methods to extract the damage sensitivity 
features necessary for visualization in the terms of pattern recognition. In this regard, it will be gathered 
from subsequent discussions that, despite its conceptual simplicity and theoretical tractability, the damage 
index possess a rather intricate vectorial/tensorial character. As a result, it is best to approach the pattern 
recognition strategy from the viewpoints of both intensity (magnitude) as well as directional properties. 
For instance, the magnitudes of various components of the damage tensorial index are shown to exhibit 
large and abrupt changes (spikes/peaks) in the presence of “true” defects. However, taken in isolation, this 
representation alone is not sufficient, since also vibration response changes due to environmental and/or 
operational variability (i.e., the so-called false alarm tests, which are known to pose extreme difficulties 
for many existing detection schemes) can trigger spurious patterns of this type. 

Hence, the directional property of the damage index is quite unique. The vector flow fields evaluated 
form the projection of the tensor component on a (variation - consistent C° - field of) position vectors 
perturbations will always be pointing in the direction of increased dissipation. In other words, vectors 
perturbations will point in the direction of decreased total stored energy. Consequently, for true defect 
case (irrespective of the underlying physical mechanism leading to the defect/deterioration) the discrete 
set of arrows representing these vector fields on the boundaries of the regions enclosing the defects will 
be directed outwards. This will persist in a consistent manner, irrespective of the mechanical response 
signature being interrogated, i.e., any vibration mode or any static load testing. On the contrary, this 
position will be completely or partially lost in case of false alarm test, e.g., all vectors will not be 
reversed, for extreme case of increase stiffness in a small localized region in the structure, to point 
inwardly. 

3.1.2 Governing Equations for Flexural Vibration of Shear-Flexible Plate. — Consider the case of 
mechanical response of a homogenous, isotropic plate (with no thermal effects, etc.) with no defects 
subjected to free flexural vibrations. Adopting the well-known approach of treating free vibration, we 
consider for one typical mode; i.e., the n th mode with frequency 00 = 00 ^ and mode shapes ( w , \|/ 1 , \|/ 2 ) 

= (w n , \| J ln , \| ; 2n )• Here w, \|/ 1 , \|/ 2 represents the displacement in z-direction and the rotations in the x 
and y- directions respectively. Note that, for convenience, we will drop the subscript V in all the 
subsequent derivations. A shear - flexible theory developed by Mindlin/Reissener(ref. 77) is utilized for 
plate flexure. This shear - flexible theory has both bending and transverse shear effects as well as lateral 
(linear) and rotary inertia accounted for. In addition, the following notation is introduced. 

The kinetic quantities, namely, the curvatures and the transverse shear may be expressed as, 


Curvatures: 


K * 



K r 
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Transverse Shears: 


dw 

T * = & +Vl ’ 
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dw 

dy 
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The static quantities, i.e., the shear resultants per unit width which taken from the moments are 
expressed as, 

Moments: 


hi 2 

M x= M U = f °xx zdz 
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h/2 
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M xy= M \2= M 2\ = j V xy zdz 
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The constitutive relationships for the plate isotropy and with the assumptions of linear-elastic 
response take the following form, 


Mu ~ D h 
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Where 
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Finally in what follows are the well-known conservation expressions in physical space. That is, the 
balance of the linear and angular momentum in the form of the first order differential equations is given 
as follows: 
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dM\ i dM \ 9 o 

— H — +p/co 2 \]/ 1 =0 (3.1.1 

dx dy 

dM 19 dM 22 9 

— (3.1- 

dx dy 

^L + ^ +P to 2 w=o (3.1.: 

dx dy 

The following sections present a step-by-step procedure from which we can show that the equations 
(3.1.1), (3.1.2), and (3.1.3) are the basic equations for the Defect Energy Parameter. 

3.1.3 The Detection Parameter Specialization . — The total strain energy, W, for thin plate can be 
defined as 


W(w,\\f u \\f 2 )~ — M n —— + M 22 + M \2 + ~ — +01 ^ + V fl + 02 ^ + ¥2 

2 ox dy y dy ox J y ox y y dy 

(3.1.4) 


2 11 0x zz av 0y dx J ~\dx ri J \dy 


After substituting the expression for M n , M 22 , M 12 , 0i and Q 2 in equation (3.1.4) the following 
expression for W is obtained 

W(w,y i,x|/ 2 ) = 

-GM 2 ^! 2 +2gM 2 ^ +GM 2 ^2f— 1 + — GM 2 f— 1 + -Z)A^l1 

2 2 y dy J 2 y dy ) 4 y dy y 




+7 V 


_ 2 D '’ V 


+ -A 


dw | 1 2 

— H — Ghk 

dw 


y dx ) 2 

< dx J 



if +' D 

Y5V2) 

v dx 

J 2 V dy 

A dx J 

f 

D b 


5^2 f 

V 

dx J 4 b 

1 

dx y 


2 , 

r 1 t 

The kinetic energy per unit area, T = —pV Vdz 

J 2 


1 1 1 • T , T d\y 2 d\ 1/1 dw 

where the velocity vector, V = x 2 , , — 

_ dt 3 dt dt _ 

After integrating over the thickness of the plate, T reduces to (for the special case of free vibrations; 
i.e., over one natural period of vibration in the typical mode number “ n ”) 


T= \ p( * n Yi h3 (^ + ^i) +hw2 
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The first term in equation (3.1.6) represents the rotary inertia due to the rotations ijq and \|/ 2 along 

the x andy-axes, and the second term is the translation/transverse inertia due to the displacement w along 
the z-axis. 

The defect energy parameter for detection can then be formed as two components of “forces”, i.e., 


F X = F \ =\{Pnl + P 2 \m)ds 

S 

F y=P 2 = J (P\2 l + P 2 2 m )ds 


> 


(3.1.7) 


Where / and m are the direction cosines of the unit normal to “ds'\ and 5 is the length of the perimeter 
curve surrounding an area of the plate middle surface, the terms P n , Pn, Pn, and P 22 are given as, 


P n =(W-T) 


Mi 


ll ' 


d \ |/ 1 
dx 


+ M 


12 


d\\f 2 


Pll=~ 


Mi 


12 


^12 =■ 


Mi 


ll 


d \ |/ 1 
dx 
d \ |/ 1 


+ M> 


22 ' 


dy 2 


+ 02 

dx dx 


_ dw 

+ Qi — 

dx dx 

dw 




,, dy 2 _ dw 

+ . V/ 1 2 h 0| 

^ dy dy dy 

P 22 =(W-T)-(m u ^ + M 22 ^ + Q 1 ^ 
V dy dy dy 


> 


j 


(3.1.8) 


In the absence of any “defects” in the plate structure, the associated forces in eq. (3.1.7) must be zero 
and equivalently this will lead to the following condition: 

Divergence of tensor P = 0. (3.1.9) 

The above “proof’ will serve to facilitate the understanding of the abstract, “mathematical-intriguing” 
nature of the damage parameter components. Note that, the formal derivation of the proper expression for 
such damage parameters requires the study of the variational symmetry of complex structures governed 
by coupled systems of partial differential equations. Using the language of “finite-strain” analysis, one 
then considers the inverse motion of the structure i.e., variations in material space while then present 
deformed state remain unchanged. 

According to the defect energy point of view, the basic parameters in the equations (3.1.7) to (3.1.9) 
convey the information about the damage complete and clear. One can consider changes in response due 
to the defect energy equations to give changes in natural frequency or mode shapes by making some 
approximations and assumptions. Though this is useful in some cases, the ensuring detection indices will 
certainly be limited in their scope for general applications. However, if there are no approximations, all 
individual contributions such as natural frequency, deformations, stresses, strains, shear and bending 
forces, etc. interact in a very complex manner with their respective different weights and degrees of 
participation on a particular excitation to produce the final damage index. Here, the latter argument is 
adopted in more comprehensively. 

3.1.4 Discussion and Generalization . — With above in mind, one can obtain the necessary 
expressions required for the analysis of other structural vibration modes; e.g., the membrane/inplane 
counterpart. In fact, a more effective approach is to consider the combined flexure/membrane formulation 
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as in the case of general, spatially curved shells. This indeed was the case considered in the 
implementation completed for the purpose of the feasibility study for the damage detection scheme. 

In particular, the details of the formulation utilized here as a basis for the numerical simulations is 
given in (ref. 80); i.e., for general mixed-type, shell finite elements. In the sequel, only a very brief outline 
of the schemes will be given. As to the simulation of damage scenarios, the simplest types of defects have 
been considered; i.e., in the form of “elastic” material stiffness degradation. Note that this approach has 
also been the most popular approach in the existing literature on damage detection. Further generalization 
to more complex types of dissipative damage/defect phenomena should be straightforward, both 
conceptually as well as mathematically. 

3.1.5 Background on Simulations. — Theoretical development of the damage index is based on the 
shear flexible theory of Mindlin/Reissner (ref. 77). This theory includes rotary inertia and shear in the 
same manner as Timoshenko’s (ref. 81) one-dimensional theory of bars. 

Depending upon the treatment of the effect of transverse shear deformation, the existing plate bending 
elements may be divided into two groups: one based on Kirchhoff Plate Theory and the other on Mindlin 
Plate Theory. In the formulation according to Kirchhoff, finite elements derived from the principle of 
virtual work or the principle of minimum total potential energy must satisfy the C 1 compatibility across 
element boundaries. On the other hand, the Mindlin plate formulation, that includes the effect of 
transverse shear deformation, C 1 compatibility is not required, even an element derived from the principle 
of minimum total potential energy. In addition, a Mindlin plate bending element can be easily extended to 
a degenerate type shell element, with curved geometry. But unfortunately Mindlin plate bending elements 
have a tendency to lock as the thickness of the plate becomes very small. Therefore, in formulating a 
Mindlin plate bending element, special care must be taken to eliminate locking. A similar locking effect 
has been discussed by Hinton, et al. (ref. 82). 

The development of suitable finite element models for linear and nonlinear analysis of plates and 
shells has always presented a major challenge, due to the many theoretical intricacies involved. To 
overcome such intricacies, Saleeb, et al.(ref. 79) developed a simple, shear flexible, quadrilateral plate 
element, designated as HMPL5. To predict the capabilities of the damage index, this quadrilateral plate 
element has been used for the present computer simulations. This element has five nodes, with three 
displacements and two rotations at each node. The interior fifth node is located at the geometrical centroid 
of the element. As Mindlin plate theory accounts for transverse shear in addition to bending deformation, 
this shear flexible quadrilateral plate element is valid for thin as well as thick plates. 

3.1.6 Outlines for Shell Modeling in Numerical Simulations. — Based on a modified Hellinger- 
Reissner variational principle, where both displacement and strain fields are assumed independent, 

Saleeb, et al.(refs. 78, 79, and 80) developed an effective 5-node shell element, designated as HMSH5. 

The HMSH5 element is primarily used for damage numerical simulations. Five degrees of freedom are 
defined at each nodal point, that is, three translations ( u , v, w) along the Cartesian global axes and two 
rotations ( Qj , 0 2 ) about mutually perpendicular lamina coordinates. In total the HMSH5 element has 25 
DOF. In finite element descritizations using element type HMSH5, the displacements are interpolated in 
terms of nodal degrees of freedom and can be written as 


u = Nq (3.1.10) 

where N is the interpolation matrix and q is the vector of nodal displacements of the element where 

nr 


q 


u l ,v l , ,M 5 ,v 5 ,w 5 ,ej J \e^ 


(5) fl(5) 


(3.1.11) 


•• 

in which the superscript T denotes transpose of a matrix. Consequently, the acceleration field, u , within 
the element is interpolated in terms of nodal acceleration as 
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(3.1.12) 


u =N q 

while strain 8 are approximated in terms of a strain parameter (3, as 


e = /?p (3.1.13) 

where (3 has 19 terms, among which seven terms that belong to the ‘membrane’ lamina strains (constant 
through thickness) and 12 terms are included in the combined bending part. In the combined bending 
portion, five terms are the transverse shear strain components. P is a (5x19) strain interpolation matrix for 
element lamina strains. In general, the entries in P are polynomial functions of natural coordinates. By 
utilizing equations (3.1.10) and (3.1.13), the Hellinger-Reissner functional n R can be written as 


n* =q T Mq-±tfHp + tfGq-W exl (3.1.14) 

Where 

H = ^P T DPdv, G = ^P T DB l dv, and M = ^ N T pNdv 

V V V 

In the above, D is the elastic stiffness matrix, B l is local strain displacement matrix and W ext is the 

external work done. Invoking the stationary behavior of equation (3.1.14) with respect to strain yields (3 in 
terms of q 


(3 = H~^Gq (3.1.15) 

which is used to eliminate the strain parameters on the element level. Finally substituting equation 
(3.1.15) into (3.1.14) results in 

n R =q T M'q+ } -q T Kq-W ext (3.1.16) 

where the stiffness matrix for the hybrid/mixed element is given by 

K = G T H~ l G (3.1.17) 

Once the stiffness and mass matrices are assembled for the entire shell, the equilibrium equation can 
be obtained as 


M q + Kq = 0 (3.1.18) 

which leads to the standard generalized eigen-problem to be solved for co and the corresponding eigen- 
vector, (|>; e.g, subspace iteration (Bathe [ ]) 


(at — <o 2 A/')<|> = 0 


(3.1.19) 
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3.1. 7 Damage Scenario . — When damage develops, effective properties of the material in the 
damaged area of the structure are changed. These changes in the material properties may be characterized 
by modification of the stiffness as well as the damping coefficient which pertains to the vibration 
behavior. In the present study, a 2-D plate was chosen for damage detection using the defect energy 
parameter since detailed investigations of other beam and frame like structures has been already 
conducted. 

For each plate, several damage locations are imposed. These can be a single or multiple situations that 
may be imposed at any span of the plate, interior or exterior, close to or away from prescribed boundary 
conditions. For each location, different amounts of damage severity may be imposed. Damage severity is 
represented by the percentage reduction of bending rigidity, El, as well as shear rigidity GA. In this study, 
the percentage of damage is defined as that specific percentage reduction of Young’s Modulus (E) at the 
specified element. For example, “ 1 percent damage at element 120” means a 1 percent reduction of 
Young’s Modulus (E), i.e., 0.99E is used in element 120 in the finite element model. Other material 
properties such as cross-section, mass and moment of inertia remain intact. Along with the percentage 
reduction in the Young’s Modulus, the damage also simulated as a line crack by creating multiple nodes 
at the same geometric location and changing the nodal connectivity, and also by eroding an element to 
create a tiny hole in the structure. 


3.2 Wavelet Theory 

3.2.1 Introduction . — The concept of wavelets has been discussed in the literature since the early 
1800’s. However, only in recent years has significant progress been made in the applications of the 
wavelets to signal processing. Historically mathematics, physics, electrical engineering, and other applied 
sciences have contributed to the development of the wavelet theory. For example, significant practical 
applications of wavelets have been found in signal and image processing, wireless communications, and 
control system analysis, which are in the electrical engineering domain. The pioneering work of 
Daubechies in the early 1980s has shown effective use of wavelets in so-called interdisciplinary research. 
Other fields that are making use of wavelets include astronomy, acoustics, nuclear engineering, signal and 
image processing, music, magnetic resonance imaging, speech discrimination, optics, earthquake- 
prediction, radar, human vision, and pure mathematics applications such as solving partial differential 
equations. 

If we consider signal processing, wavelets are used for transforming a signal into a form that it is 
more useful to the application. For example, if we are interested in reducing the noise in a signal, the best 
representation is the one in which the signal and noise are easily separated. The two most common 
representations for a one-dimensional signal are the time signal, and its Fourier transform. Unfortunately, 
it is not easy to extract frequency information from the time signal and vice versa. Different signal 
representations depend on the choices of the time-frequency resolution. The choice of the proper signal 
processing technique is based on the signal that needs to be analyzed. Short-time Fourier transform 
(STFT) and the continuous wavelet transform (CWT) are used to represent a signal in both the time and 
frequency domains. 

3.2.2 Continuous Wavelet Transforms . — The wavelet transform maps a function, s(t), into a two- 
dimensional domain and is denoted by 


W s {a,b) = 



a J 


dt 


+00 

J s(t)Kib { { ) dt 


—00 


(3.2.1) 
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where h(t ) is in general called the mother wavelet, and the daughter wavelets, are given by: 


uo 



ft-b\ 

a J 


(3.2.2) 


Equation (3.2.1) is also known as the forward transform or analysis; h ab ( t ) is a set of basis 
functions/daughter wavelets obtained from the mother wavelet h(t) by compression or dilation using the 
scaling parameter a and translation using the shift parameter b. The scaling parameter a is positive and 
varies from 0 to oo. For a < 1, the transform performs compression of the signal, and for a > 1, the 
transform performs dilation of the signal. The reason for choosing the factor l/ a [a in the above equation 
is to keep the energy of the daughter wavelets constant. Without this normalization factor, for different a 
values the wavelets dilate or compress and their total energy changes. 

The inverse wavelet transform is given by: 


+00 00 


s(t) = ~ J \w s (a,b)h 


t -b 


-oo 0 


y a J 


da „ 
—~db 

a 2 


(3.2.3) 


in which the constant c is defined as 


c = 


i 


—00 


2 

—dco<co 
co 



(3.2.4) 


Equation (3.2.3) is also referred to as the reconstruction formula, inverse transform, or synthesis, and 
equation (3.2.4) is known as the admissible condition. 

The wavelet theory as applied in the present research, is utilized for the reconstruction of the damage 
“signal” from the noisy data. In this context, it is necessary to treat the damage signal as a discrete signal. 
Therefore, discrete wavelet transforms are used to complete the study. The following section describes the 
theory behind discrete wavelet transforms in detail. 

3.2.3 Discrete Wavelet Transforms. — A discrete signal is a function of time with values occurring at 

discrete instants. Generally we shall express a discrete signal in the form / = (/ , ^ j fi * , fn ) > 

where N is a positive even integer which we shall refer to as the length of f. The values of f are generated 

by sampling a continuous signal g at uniformly spaced time intervals, i.e., t-t x ,t 2 , , t N . The values of 

/ are 


/l ~ g(h X fi - g(t 2 )’ >/ n ~ S(. { n) (3.2.5) 

All wavelet transforms decompose a discrete signal into two subsignals. Each of these subsignals is 
half the length of the original signal. One subsignal is referred to as the average or trend while other 

subsignal is referred to as the difference or fluctuation. The first trend subsignal, a = (a x ,a 2 , , a N/2 ) > 

for the signal /is computed by taking the average of the values in f. The first fluctuation of the signal f 

which is denoted by d l =(d l ,d 2 , , d N / 2 ) , is computed by taking difference of the values in f. 

To further clarify the concept of the discrete wavelet transform (DWT), the Haar wavelet, which is 
the simplest type of wavelet, is taken as an example and will provide a good foundation to understand the 
more sophisticated wavelet transforms. 
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3.2.3. 1 Haar Wavelet Transform: The Haar transform of the signal / = (/i 9/2 5/3? >/#) w ill 

decompose the signal into averages, a , and differences, d. The decomposition can be done at multiple 
levels. The first level is the mapping H 7, defined by 




f — K# 1 1^ 1 ) 


Averages: a x = (aj , a 2 , • 


•> a N!2 ) 


aj =■ 
a 2 = 


(/1+/2) _ yr( (/l +/ 2 ) 

V2 

(/ 3 +./ 4 )_^ 


2 

'(/ 3 +/ 4 ) 


V2 


v 




a N/2 


(/iV-1 +/jv) 


>/2 




(/am + /v) 


Fluctuations: d x ={d\,d 2 , ,d N / 2 ) 

f/j =- 


'A 


d 2 = 


(/ 3 -/ 4 ) 




s 

(Ox-fi) 
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\h ~U) 






7V/2 


C/:V ~1 ~ /Af ) 


=v^ 


C/iv -1 - /v) 


Inverse Haar transform from level- 1 decomposition: 


(3.2.6) 


(3.2.7) 


(3.2.8) 


(</' )<:/')- 


— £/] (^2 + ^2 — ^2 


a N/2 +c ^N/2 a NI2~^NI2 


V2 ’ V2 ’ V2 ’ V2 


V2 


V2 


= / (3.2.9) 


The second level Haar transform decomposes the trend subsignal a . The second level mapping H 2 , 
defined by 


f^(a 2 \d 2 \d x ) (3.2.10) 

And the N level mapping H N , defined by 

f ->(a N | d N | \d 2 \d l ) (3.2.11) 

3.2.3. 2 Haar wavelets and scaling signals : In this section we further discuss the Haar wavelets and 
the scaling signals, and defining the averages and fluctuations as an inner product of these signals. The 1- 
level Haar wavelets are defined as 
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Consider the inner product of / and w\ , 



Therefore, we can summarize that, 1 -level differences are defined as, 

d l =( f W 1 ) 

u m ’ rr m j 

in which the 1 -level scaling signals are defined as 



In addition, consider the inner product of f and v\ as, 



Therefore, we can summarize that, 1 -level averages, are given as 


(3.2.12) 


(3.2.13) 


(3.2.14) 


(3.2.15) 


(3.2.16) 
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(3.2.17) 


a 1 =lfV X \ 

’m / 

Thus in general, 

Averages, a% =(^f 9 V^ and Differences, d% = (V, ^ (3.2.18) 

and the inverse Haar transform as an inner product is defined as, 

f = A k +D k + + D 2 +D 1 (3.2.19) 

wher -e, A k =(f,V k )v k and D k =[f,W k )w k 

Equation (3.2.19) shows how a discrete signal is synthesized by beginning with a low resolution 
signal and successively adding on “details” to create a high resolution signal, ending with a complete 
synthesis of the signal at the finest resolution. This process is known as “multiresolution analysis 
(MRA)”. 

The typical application of the Haar transforms is compression of a signal and removing noise from a 
signal. In this report, wavelets are used to reduce the noise intensity in a signal. Therefore, we will focus 
this particular application and the following sub-section describes on the process of noise removal using 
the Haar wavelet. 

3. 2. 3. 3. Denoising of a Signal Using Haar Transforms : A basic signal model gives the measured 
signal as actual signal plus its noise, i.e., 

/= s + n 

where / = Measured signal 

5 = Actual signal 

n = Noise 

Note the term noise is referred to as the disturbance in the measured signal. To remove or to reduce 
the intensity of the noise in a signal, the following procedure is used: 

(1) Take the wavelet transform of signal f 

(2) Choose a threshold T, so that the noisy signal’s transform values have magnitudes which lie 
below T. 

(3) Set values of the transform to zero which are less than the threshold T. 

(4) Take the inverse transform to get the denoised signal. 

In this research, Daub4 wavelet from Daubechies family is utilized. The Daubechies wavelet 
transforms are defined in the same way as the Haar transform by computing averages and differences via 
scalar products with scaling signals and wavelets. The only difference is how the scaling signals and 
wavelets are defined. 

3. 2.3 A TheDaub4 Wavelets and Scaling Signals: In this section, details of the Daub4 wavelet are 
presented. First, the scaling numbers a 1? a 2 , a 3 a 4 are defined by 


1 + V3 _3 + V3 

W2 ’“ 2 ~ W2 1 


1-V3 _3-V3 

W2 ’ a “~ W2 


(3.2.20) 
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Using these scaling numbers, the 1 -level Daub4 scaling signals are expressed as, 


v \ = (a 1 ,a 2 ,a 3 ,a 4 ,0, 0, ,0) 

V\ = (0, 0, oq , a 2 , a 3 , a 4 , 0, 0, ,0) 

Vl =(0,0,0,0,a 1 ,a 2 ,a 3 ,a 4 ,0,0, ,0) 






V X N =(0,0, ,0,ai,a 2 ,a 3 ,a 4 ) 

T _1 

V X N = (a 3 ,a 4 ,0, 0, ,0,a 1? a 2 ) 

T J 


Let us define, 


^=( 1 , 0 , , 0 ) 

F 2 °= (0,1,0, ,0) 


f 2 °=(0,0, ,0,1) 

Using the equation (3.2.22), the first level Daub4 scaling signals satisfy 

V m — (Xj V2 m -\ + ^4^2w+2 

Similarly, the second level Daub4 scaling signals are defined by 

V m = a l V 2m-l + a 2 V 2m + a 2 V 2m+\ + a 4^2m+2 

All other levels of Daub4 scaling signals are defined in a similar fashion. 
Let the wavelet numbers Pi , P2 , P3 , P4 be defined by 


Pl=- 


V3 _ _V3 

’ P2 “ 


4V2 


4V2 


3 P,=4#,f> 4 = 


4V2 


-1-V3 

4V2 


Notice that the wavelet numbers are related to the scaling numbers by the equations: p] = a 4 
P 2 = -a 3 , P 3 = a 2 , P 4 = -aj . Using these wavelet numbers, the 1 -level Daub4 wavelets 

W\ , w \ , , Wtf / 2 are defined by 
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(3.2.21) 


(3.2.22) 


(3.2.23) 


(3.2.24) 


(3.2.25) 



fFi 1 =(pi,p 2 ,P 3 ,P4,0,0, ,0) 

= (0, 0, Pj , p 2 , p 3 , p 4 , 0, 0, ,0) 

W 3 = ( 0 , 0 , 0 , 0 , Pj , p 2 > P 3 5 P 4 5 o? 0 , , 0 ) 




K =( 0,0 , 

T _1 

^=(p 3 ,p 4 ,0,0, 


? 0 , Pi , P 2 » P 3 ’ P4 ) 


■APl,p2) 


2 


J 


The 1 -level Daub4 wavelets satisfy 


K =Pl^Vl + P 2 V 2 °m + P 3 ^+1 + P 4 4 )+ 2 
Similarly, the 2-level Daub4 wavelets are defined by 

K 2 , =Pl^-l +P 2 ^2 W + WL + 1 +p4^ + 2 


(3.2.26) 


(3.2.27) 


(3.2.28) 


All other levels of Daub4 wavelets are defined in a similar fashion. 

3.2.4 Application of Wavelets in De-noising Data . — As mentioned, wavelets are being applied in a 
variety of applications. A few of the applications are in computer and human vision, FBI fingerprint 
compression, denoising of noisy data, detecting self-similar behavior in a time-series, and musical tones. 
In recent years researchers started using wavelet theory as a technique for damage detection. But in this 
study, it is used as tool for denoising the measurement data. 

In this study, wavelets are utilized in the signal feature reconstruction from the noisy data. When we 
decompose a signal using wavelets, we use filters that serve as averaging filters and others that produce 
details. Some of the resulting wavelet coefficients correspond to details in the signal. If the details are 
small, they might be omitted without affecting the main features of the signal. The idea of thresholding is 
applied to set to zero all coefficients that are less than a particular threshold. These coefficients are used 
in an inverse wavelet transformation to reconstruct the signal. 

There are two types of thresholding used in denoising, i.e., hard thresholding and soft thresholding, 
figures 3.1 and 3.2 shows the hard and the soft thresholdings, respectively. Hard thresholding is not a 
continuous function, so it will exaggerate the gaps between significant and insignificant values. Soft 
thresholding is a simple modification of hard thresholding such that is becomes a continuous function 
which as a result, it does not exaggerate the gap between significant and insignificant values. In general, 
the choice of thresholding depends on the type of application. 


(1) Hard Thresholding 


H{x) = 


x if | x | > T 
0 if|x|<r 


(3.2.29) 
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T 

0 

- T 


Figure 3.1. — Hard thresholding. 



T 

0 

- T 


Figure 3.2. — Soft thresholding. 


(2) Soft Thresholding 


S{x) = x if | x | > T 

2 x-T if 772 < x<T 
T+ 2x if -T<x < -772 
0 if \x\<T 


A 


7 


V 


J 


(3.2.30) 


Let us consider a noisy signal, which is taken from the wavelet toolbox demos and it is called as noisy 
blocks, as an example to see the effect of denoising by using a fifth level Daub4 wavelet found in the 
Wavelet Toolbox in Matlab [ ]. The Daub4 wavelet is denoted by “db2” in wavelet toolbox. Here, 
different values of soft thresholding, T 7 , are used. Figures 3.3 to 3.6 show the effect of the various values 
of the threshold, T. As expected, we can observe more reduction of the noise in the signal by increasing 
the value of T. The values of the T need to be “optimized” in order to retain the significant values of the 
original signal. Figure 3.7 shows the details/differences of the noised signal with a threshold, T= 4.219. 
The dotted line in each level of details (d\ to d 5 ) represents the value of each threshold in which the 
denoising coefficients below that threshold value are set to zero. 
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Figure 3.3. — Signal with the noise. 



Figure 3.4. — Denoised signal using threshold = 0.5. 



Figure 3.5. — Denoised signal using threshold = 1. 



Figure 3.6. — Denoised signal using threshold = 4.219. 
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Figure 3.7. — Details/differences of noised 
signal with threshold = 4.219. 


3.2.5 Specific Application of Wavelets in Global Damage Detection . — The signal considered in the 
above is an example discrete signal in which a scalar quantity is a function of time (as in most typical 
applications of wavelets). In this present research, the damage parameter values (a scalar) are interpreted 
as a discrete “signal”, which is not a function of time, but instead is a function of position. That is, the 
scalar damage parameter is calculated at various points in a structure, which for the case of a flat plate 
have a specific position (x,y) value associated with it. Thus it is necessary to convert the “two- 
dimensional” (spatial) array of damage parameter values to a one-dimensional array of values. To achieve 
this, the values of the damage parameter are processed in a sequential manner over the area of the plate. 
Figure 3.8(a) shows the typical sequence in which the damage parameter values are read. It is very 
important to note that when forming the damage parameter signal, the sequence in which the values are 
read is of prime consideration. 
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Figure 3.8. — A plate with a mesh size of 13 by 16, wherein three elements are damaged 
(a) Sequence of damage parameter consideration is “left to right” (b) Sequence 
of damage parameter consideration is “bottom to top.” 




Figure 3.9. — Damage parameter signal (a) Sequence from “left to right” 
(b) Sequence from “bottom to top.” 


Figure 3.8 shows a plate with a mesh size of 13 by 16. In the total of 208 elements, three elements are 
damaged at a location, i.e., treating as single damage in the plate. When the damage parameter is 
generated, damaged elements will have the maximum values. When we form a discrete signal by 
considering the damage parameter values from “left to right” as in figure 3.8(a), then the signal will have 
peak at single location as in figure 3.9(a), since three damaged elements are next to each other in this 
sequence. But, if we consider from “bottom to top” as in figure 3.8(b), then the signal will have peaks at 
three locations as in figure 3.9(b), since three damaged elements are not next to each other in this 
sequence. In both cases, there is single damage. Therefore, in the signal feature reconstruction, care has to 
be taken on the sequence of damage parameter when forming a discrete signal. 

In the following sections of this report, a simulation case study is performed pertaining to signal 
feature reconstruction using the discrete wavelet transforms outlined above. The wavelet toolbox in 
Matlab is used for all simulations and the results of this case study are presented in chapter V. 
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Chapter IV 
Computational Tools 


4.1 General 

In the general engineering community, the use of computational simulation software has become an 
obvious solution for many of the challenging engineering problems. In engineering, the design/analysis 
phase requires an iterative process in order to obtain the optimal design. Computational simulation 
software can handle this iterative process very efficiently, so an engineer can concentrate on other major 
issues of the work. In this context, many software tools have been developed for engineering applications. 
Software tools such as ABAQUS, ANSYS, LS-Dyna, etc. are examples of widely used general purpose 
finite element based structural analysis programs. Other tools such as Matlab, Mathematica, etc. are 
useful for such topics as signal processing, optimization, symbolic derivation, etc. All of these tools serve 
to increase the productivity/efficiency of the engineer. 

With the above points in mind, an efficient computational tool to assist in structural damage detection 
is developed in this research. A computational tool, which will be referred to as the D-Tector, previously 
developed in (refs. 67, 68, and 85), is used for the global (structural) damage detection. With D-Tector 
serving as the central software component, other software packages such as ABAQUS, Visual C++, 
Matlab, and Wavelet Toolbox are combined in a computationally efficient manner to allow for the 
systematic study to be presented in subsequent chapters. Figure 4.1 describes how these tools were 
utilized. The following sections present a detailed explanation of the implementation of these tools in this 
research. 


4.2 D-Tector for the Global Damage Detection 

4.2.1 Introduction . — Observing the necessary usefulness of having a capability for structural 
diagnosis and prognosis, this present research study focuses on using a finite element based algorithm for 
the damage detection and localization of a structural member. Primarily, this study concentrates on the 
development of a robust scheme for the global damage detection which utilizes as a central component in 
the D-Tector software. 

The main purpose of D-Tector is the detection, location and determination of the approximate size of 
defect/defects in a structural member, using the displacements measured during the experiments. 
Considering the use of the latest experimental techniques, viz. the Electronic Speckle Pattern 
Interferometry (ESPI) method of producing interferograms, D-Tector contains a feature to convert the 
experimental file written in binary format to ASCII format and ultimately to the input file format 
necessary for the damage detection core. This computational tool contains a series of programs written in 
Fortran and C and a graphical interface, called Eview, for post processing. Eview is based upon various C 
Graphics libraries and is used visually interpret/extract the damage. 

4.2.2 Structure of the D-Tector . — A program’s efficiency is directly related to the structure of any 
computational tool. The D-Tector tool has three major sections: Preliminary Processing, Damage 
Detection Core, and Post Processing. Each section contains its own subroutines with a special purpose. 
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Visual C++ 
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for viewing results 
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RESULTS 


Structural Diagnosis and Prognosis 


Figure 4.1. — Computational tools for global damage detection. 


D - Tector 



Figure 4.2. — Structure of the D-Tector. 


The developed NDE scheme for global damage detection is contained in the Damage Detection Core 
section. Therefore, this section will handle the necessary computations for the detection. As written, the 
damage detection core can understand only specific formats of the input, but the experimental files may 
be in a binary format which is not compatible. The process to modify the experimental files into a format, 
which is compatible to the damage detection core, is undertaken by the Preliminary Processing section. 
The format of this input file has been designed in such a way that the user can generate it manually, if the 
frequency of vibrations and the displacements of the structure at specific locations are known. In this case 
the so-called preliminary processing is not needed. However, since many experimental techniques today 
are conducted using digital computers for data acquisition, and which may be in binary format, the pre- 
processing capabilities are required. Similarly, the post-processing section is for viewing the files 
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generated by the damage detection core for the so-called “feature extraction.” By viewing the files with 
Eview, the user can detect and locate and measure the approximate size of the damage. Eview also has the 
provision for plotting the mode shapes of the structure. 

4.2.3 Algorithm of the D-Tector . — The initial stage of the D-Tector is the preliminary processing. It 
contains four applications viz. Lazout2surfit.exe, Surfit.exe, Pix2ev.exe, and Pix2evNonorm.exe. The 
input file for the damage detection core will be generated from the experimental file in this section. The 
figure 4.3 illustrates the series of the applications utilized. 

A typical laser holographic testing experimental setup that measures the modal response, which is 
used for the global damage detection, generates a binary output file with .out extension. To retrieve the 
holographic data, i.e., to convert the binary file to ASCII format, the Lazout2surfit.exe program is 
utilized. This program takes the experimental output file as input and coverts into the surfit file, which is 
in ASCII format. The retrieved holographic data in the surfit file appears as an array of pixel values. Also, 
reduction of the data may be required if large volumes of data is produced for the efficiency of the 
analysis. The surfit.exe program is used to translate and truncate the data, i.e., for ignoring particular 
regions of data. This code fits the translated data and samples the data at a resolution that the damage 
detection core can use. The translated data from the surfit.exe is stored in the .xyz file. The .xyz file 
consists of number of columns, which indicates the location and the corresponding intensity at that 
location of the plate according to the dimensions of the plate. The first column indicates the x-coordinates, 
second column indicates the y-coordinates and the third column indicates the intensity/value at that 
particular location. The sequence continues until the end of the plate. 



Figure 4.3. — Flow chart for preliminary process. 
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The damage detection core accepts a specific format of the input file. Pix2ev.exe generates the 
damage detection core input decks with normalized displacements by taking the ud.xyz (undamaged) file 
and the d.xyz (damaged) file as input. Also, there are few variables that need to be input in addition to the 
damaged and undamaged files, such as the frequencies, values for the mesh size, and the thickness of the 
plate. There is the option to put this data in an “answer” file so that the program may be run in a “batch 
mode” as opposed to interactively to decrease the burden on the user. All of the input variables are in the 
answer file, so the user only has to change the damaged and undamaged frequencies, values for mesh size, 
and the plate thickness for each analysis. Pix2ev then creates a data input (.din) file, which will be used as 
the input file to the damage detection core. Figure 4.4 illustrates the major segments of the typical input 
file. The data input file consists of the degrees of freedom, nodal and element definitions, boundary 
conditions, frequencies of vibration of the undamaged and damaged plates, and the undamaged and 
damaged displacements. 

The displacements are normalized when pix2ev.exe is adopted for developing the data input file. 
Pix2evNonorm.exe is also develops the data input file (.din) but without normalizing the undamaged and 
damaged displacements. If you want to consider the actual values instead of normalized values in the 
detection this application is needed. 

The data input file generated by the preliminary processing will now be in a format that is usable to the 
damage detection core. At this point, the data will be sent to the damage detection core section where the 
damage computations are performed. The damage detection core section generates the files in such a way 
that the user can view in Eview the location and approximate size of the damage. This section contains 
three main applications for the D-Tector, viz., asg_dam.exe, asg_get_dam_vector.exe, and evdelta.exe. 

Figure 4.5 emphasizes the sequential list of applications in the Damage Detection Core section. The 
asg_dam.exe program is the heart of the D-Tector. The major formulation for the damage detection is 
coded in this application. It has the ability to perform the static and dynamic analyses and the damage 
parameter computations. Specifically, this code is finite element based. Since this study focuses on plate- 
type structures, a 3-D shell element with five degrees of freedom (3 translations and 2 rotations) is used 
primarily. For generality, the code also has ability to work with four types of elements viz. truss, beam, 
curved beam, and shell elements. Asg_dam.exe produces the output (.dot) file. This data output file is 
very large and very significant in the detection process. This file contains nodal and element definitions, 
displacements, stress calculations, damage parameter calculations, load case data, and total system data 
i.e., information about the number of equations, the number of matrix elements and the bandwidth. By 
simply reading large “raw” data output file, it is very difficult to extract the significant damage features. 
Thus, the asg_get_dam_vector.exe program is used to extract the specific information required to 
graphically view the results. It generates five files viz. raw damage parameter (.rdp), smooth damage 
parameter (.sdp), damage vector (.dav), undamaged mode shape (.uds), damage mode shape (.dds). The 
.rdp file contains the raw nodal damage parameter vector definitions and element data damage 
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Figure 4.4. — Major segments in data input file. 
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Figure 4.5. — Flow chart for damage detection core. 

parameter definitions where as the .sdp file contains the smoothened nodal damage parameter vector 
definitions and element data damage parameter definitions. The nodal damage parameter values are per- 
element values that have been area-weighted at the nodes. The .dav file contains the damage force 
vectors. These are calculated based on the energy vectors in the structural member due to the vibrations. 
The undamaged and the damaged mode shape files contain the displacements extracted from the data 
output file. All of these files can be viewed in the post-processing section. 

The damage detection core has the option to deal with either dynamic or static analysis. Thus, the 
asg_dam.exe program also has a provision to accept different formats of the data input file. We can feed 
the displacements which are generated experimentally or analytically. Analytical data can be generated by 
using any finite element analysis software. If we provide the displacements to the damage detection core 
then D-Tector uses the series of programs listed in the figure 4.5. In this case, the analysis is treated as 
static. The damage detection core can generate the dynamic response for the purpose of simulations. For 
the dynamic case the input file is different from the static case. When testing analytically we must provide 
the number of vibration modes with damaged element/elements definitions since the displacements are 
unknown. The code will generate the displacements according to the out-of-plane or in-plane action of the 
analysis. If performing purely a static analysis the user should provide the applied force definitions on the 
nodes. The damaged elements can be defined by the degradation of the Modulus of Elasticity . The 
damage calculations are performed and written in the data output file. In the analytical case to acquire the 
mode shapes a separate program, asg_get_mode_vector.exe, was developed. It generates undamaged and 
damaged mode shape files with .mdv extension. 
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The displacements of the undamaged and damaged plates are different in value. This difference is 
because of the presence of the damage. Therefore, if we evaluate the difference between both 
displacements, we identify that the delta mode shape is completely different than the damaged and 
undamaged mode shapes. This delta mode shape is a measure of the accuracy of the experimental 
technique. The evdelta.exe program is developed to serve this purpose. The undamaged and damaged 
mode shapes generated either by the analytical or by the experimental files will be sent to the evdelta.exe 
as input to get the difference. This file has the .dms (delta mode shape) extension. 

The final and most important section in the D-Tector is post-processing. The files generated by the D- 
Tector contains huge amount of data and identifying the significant damage values is a difficult task. A 
viewer is necessary to extract the feature of the damage in a proper pattern. In conjunction to this, Eview 
is developed. Figure 4.6 shows the Eview window. The files generated by the damage detection core can 
be viewed in Eview to locate the damage. It has developed using various C graphics libraries and has the 
feature of NURB sampling fidelity for interpolation, which allows the better viewing of the contours to 
help in localizing the damage. The mesh of the member can be viewed with nodal and element numbers 
displayed. 

The E-view has provided six different views to extract feature of the damage for specific purposes. 

By selecting the first option user can view the picture in black and white shades according to the values. 
The second option differentiates the positive and negative values. The contours are plotted in the colors, 
positive is white and negative is red, for easy differentiation. Positive values are indicated in white to 
black shade and negative values are indicated in light red to dark red shade in the picture. The third option 
uses the color spectrum to indicate the values in the picture. The first three options use the NURB sample 
fidelity and viewed as contours. Fidelity can be increased to sharpen the picture. The fourth option shows 
the magnitude per element. The fifth option is to view the values of the files in the vector field. There is a 
provision made to change the vector view threshold to hide the unwanted vectors to extract the damage. 
The vectors that are below the threshold value are erased from the picture to indicate the damage location 
more clearly. There is also a provision is made to increase the scale of the vectors. The final option is to 



Figure 4.6. — Eview window. 
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view the files in a meshed pattern using rainbow colors to indicate the values. Depending upon the type of 
file, the user can choose the appropriate contour to extract the required feature. The figures presented to 
describe the results of this study were produced by the E-View. Different types of contours produced by 
the E-View are described in the further chapters. 

The features that can be extracted from the Eview are presented from figures 4.7 to 4.12. Here, Mode- 
5 is taken as example to show the features. 
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Figure 4.7. — Damage parameter contour using black and white intensities. 
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Figure 4.8. — Mode Shape using different colors for positive and negative directions. 
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Figure 4.9. — Damage parameter contour using rainbow colors with NURB fidelity. 
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Figure 4.10. — Damage parameter contour using rainbow colors without NURB fidelity. 
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Figure 4.11 . — Damage vectors. 



Figure 4.12. — Mode shape. 


In summary, the main computational tool for detecting the damage is D-tector. Up to this point, we 
have not considered the possibility of noise in the experimental measurements. Since the primary focus of 
this research is for extracting the damage in the presence of noisy data measurements, the following 
sections will describe the utilization of other software tools necessary to process the data. 
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4.3 Visual C++ for Connectivity and to Study the Effect of Noise 


In reality, obtaining noise free measurements are very difficult. Therefore, assessment of the accuracy 
of the measurement data to noise is needed for a robust detection scheme. To simulate the actual 
laboratory environment, computer generated random numbers are added to the data sets. Here, Visual 
C++ is utilized to generate and add the noise in which the “rand( )” function in the Visual C++ library is 
used to generate the random numbers. It also used to connect the D-Tector to other tools like ABAQUS 
and Matlab. The main applications developed by the Visual C++ are GK_Noise_Damage_Parameter.exe, 
GK_InputFromAbaqus.exe, GK_GenNoise.exe, and GK_AddNoise.exe. 

GK_Noise_Damage_Parameter.exe connects the D-Tector to Matlab and loads the noised damage 
parameter values to Matlab and extracts the “denoised” values from Matlab. GK_InputFromAbaqus.exe 
connects the D-Tector to ABAQUS and will extract the mode shapes generated by the ABAQUS modal 
analysis and creates the data input file in form that is compatible to D-Tector. GK_GenNoise.exe 
generates the noise, random numbers, and writes into a file. GK_AddNoise.exe adds the noise to the data 
sets in the data input file by using already generated noise files from GK_GenNoise.exe. The source 
codes of the developed applications by Visual C++ are presented in the appendix. 


4.4 ABAQUS for Detection Scheme 

ABAQUS is a general purpose commercial Finite Element Analysis software for linear and nonlinear 
analysis. It is mainly used to perform modal analysis. Mode shapes/displacements are saved in a file from 
the ABAQUS viewer. This data is used in D-Tector or wavelet toolbox wherever required. 



Figure 4.13. — ABAQUS and Matlab in global damage detection. 
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4.5 Matlab for Detection Scheme 


The Matlab is also commercial multipurpose numerical analysis software for engineers. It is useful 
for coding and plotting mathematical functions, and to perform other mathematical 
simulations/calculations. In addition, Matlab can connect to special purpose toolboxes viz., Signal 
Processing, Image Processing, Wavelet, Curve fitting, Control Systems, Fuzzy Logic Toolboxs, etc. In 
the present research, Matlab is used to connect to the Wavelet Toolbox which will be used to “denoise” 
the measurement data. 


4.6 Wavelet Toolbox for Signal - Feature Reconstruction 

In the last two decades, extensive research has focused on the practical application of wavelet theory 
to a variety of engineering problems. According to the theoretical background described in Chapter III, 
researchers have developed number of wavelets. The most popular wavelets are assembled into one 
computational tool, called Wavelet Toolbox developed by Mathworks, Inc. This toolbox contains 15 
different wavelet families with a Graphical User Interface (GUI). It contains 1-D and 2-D wavelet 
analysis with Discrete, Continuous and Stationary Wavelet Transforms. The main applications of the 
wavelet toolbox are Denoise and Compression of signals and images. In this study, it is used for feature 
reconstruction of noisy data. The Daubechies2 (Db2), which is referred as Daub4 in the literature, wavelet 
is used from the predefined wavelets in wavelet toolbox. 



Figure 4.14. — Matlab and Wavelet Toolbox for 
feature reconstruction in global damage detection. 
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Table 4.1 summarizes the purpose of the described computational tools in this chapter. The results of 
simulations performed and the application of these computational tools are discussed in chapter V. 


TABLE 4.1.— COMPUTATION TOOLS 


Computational tool 

Purpose 
in general 

Purpose 
in this research 

D-Tector 

Global damage detection 

Global damage detection 

Visual C++ 

Programming 

Connecting D-Tector to other 
tools and to study the noise 

ABAQUS 

Finite element analysis 

Performed required simulations 

Matlab 

Coding and plotting software for 
engineers 

Used to connect Wavelet 
Toolbox 

Wavelet Toolbox 

Wavelet analysis 

Signal - feature reconstruction 
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Chapter V 

Simulations and Results for Case and Comparitive Studies 

5.1 Introduction 

Damage detection is a challenging problem that is under vigorous investigation by numerous research 
groups using variety of analytical and experimental techniques. The damage detection procedure for a 
structure depends upon the level and extent of the damage, available knowledge concerning the ambient 
dynamic environment, sophistication of available computing resources, complexity of the detection 
scheme, and depth of knowledge concerning the failure modes of the structure. 

One of the most important criteria to be fulfilled in any damage detection scheme based on a 
theoretical/computational model is the robustness of the method to real applications; i.e., using data from 
a laboratory setup or field measurements. In particular, of utmost importance for the effectiveness of the 
method is its ability to determine the location and severity of the damage when the experimental 
measurements (e.g., strains, deflections, rotations, etc.), which may contain various environmental factors 
and human errors, are provided. Simulations have been done to study the effectiveness of the proposed 
detection scheme based on the following considerations. Damage in a structure can appear as weakened 
material (e.g., in low cyclic fatigue), line cracks (e.g., high cyclic fatigue), or corrosion. If we consider the 
sensor network in an experimental setup, it can be an extensive network or coarser network. The detection 
scheme should be able to determine the location and severity of any type of failure in the structure even 
with a coarse sensor network. And also scheme should withstand some amount of noise. Along with the 
above simulations, feature reconstruction methods are considered to reduce the noise intensity in the 
measurement data. Wavelet transforms are utilized for feature reconstruction in this study. The 
background on the simulations is described in detail in the following subsections. 

5.1.1 Failure Modes! Type of Damage. — Keeping an eye on the real time environment, simulations 
are performed to complete this research. In damage detection, failure modes of a structure are one of the 
important factors. The types of common failure modes are classified as Low Cycle Fatigue, High Cycle 
Fatigue, and Corrosion. The damage in a structure can occur after a low number of cyclic loadings in the 
Low Cyclic Fatigue failure mode. This type of damage is simulated by a sudden degradation of the 
Modulus of Elasticity of element/elements in the finite element model. If the damage occurs in a structure 
after a high number of cyclic loadings then it falls under the High Cyclic Fatigue failure mode. This type 
of damage is simulated as a line crack by creating multiple nodes at the same geometric location and 
changing the nodal connectivity. Corrosion of the material is also a type of damage in the structure. This 
type of damage occurs in very old structures. When the corrosion of a structure starts, then we can 
observe tiny holes. This is simulated by eroding an element from the structure’s geometry. The finer the 
mesh, smaller the hole size is. Therefore a very fine mesh is generated to simulate a tiny hole in the 
structure. 

5.1.2 Sensor Network in Measurements. — Density of the sensor network can affect the detection 
scheme. Arrangement of the sensors for the measurements is the most important factor. The 
extensive/fine network assumes an ideal, “unlimited” ability to place and measure response at all regions 
in the structural domain. However, from a practical viewpoint, the case of coarse networks, where no 
direct measurements are made at or very near to the (“presumably unknown”) locations of the damage, 
are very crucial in judging the potential of any damage detection methodology. To emphasize this point, 
all selected coarse networks completely avoid any measurements at or very near to areas of defects. To 
study the affect of sensor network, initially, a 40 by 40 mesh is considered as the extensive sensor 
network. To simulate the coarser sensor network, 20 by 20 and 10 by 10 meshes are generated and 
equivalent measurements are taken form the 40 by 40 mesh analysis. 

5.2.3 Effect of the Noise in Measurements. — In a real time environment, obtaining ideal (noise free) 
measurements are very difficult. There is always random amount of noise present in the measured data. 
Therefore, the detection scheme should be robust enough to withstand certain level of unavoidable noise. 
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It is very difficult to identify the frequency of noise in the data. According to the literature, white noise 
and pink noise can cover most of the noise frequency ranges that are typically encountered in many 
measuring devices used (e.g., accelerometers, strain gauges, etc.). Therefore, white noise and pink noise 
are considered to study the effects of noise in the damage detection scheme. White noise contains an 
equal amount of energy per frequency band. For example, the amplitude of the sound will be same from 
400 to 500 Hz as it will be from 30100 to 30200 Hz. White noise is considered as bright noise and can be 
generated by a random number generation function. Pink noise contains equal amount of energy per 
octave (span of eight notes) band. For example, the amplitude of the sound will be same from 100 to 
200 Hz as it is from 200 to 400 Hz or 10000 to 20000 Hz. The ratio of frequency of the highest note to the 
lowest note in an octave is 2: 1 . Pink noise can be generated by filtering (averaging) the white noise with a 
frequency of 8. 

In the present study, noise is added to the displacements/mode shapes to simulate the real time 
environment. To this end, computer routines (written in Visual C++) are used to generate and add the 
noise. The following equation is used to add the noise. 


U’/(°rD’/) = U l (orD l ) 


1 + N, 



D i ~ u i 

11 

j 

U j (orDj ) )_ 


where Uf - undamaged displacements 

D” - damaged displacements 

Ni - noise level 

| D t - Ui | - Norm of difference of damaged and undamaged displacements 

5.1.4 Signal - Feature Reconstruction . — With an eye towards reducing noise intensity and their 
subsequent adverse effects on the damage detection capability, different types of filters have been 
developed. The important factor in these so-called “denoising” schemes in the damage detection is 
retaining the significant, damage-produced features, of the signals. In the last two decades there has been 
extensive research on the use of wavelet theory. One of the applications of the wavelet theory is 
denoising. In this research, one dimensional wavelet analysis is used to reduce the noise intensity. 

To study the signal feature extraction capability in the global damage detection scheme, noise is 
added to the data sets. Two types of data sets are considered; i.e., (i) raw displacements, and (ii) processed 
damage parameter values. The noise level is increased until the damage parameter contour is filled with 
the noise. The damage parameters with maximum level of noise that can cover the damage are used to 
extract the feature. The damage parameter values are processed using a one dimensional discrete wavelet 
analysis to reduce the noise intensity. Daubechies2 (db2) wavelet in the wavelet toolbox, which is in 
Matlab, is used to reconstruct the feature of the damage parameter. Visual C++ routines are used to 
interface with Matlab. 

5.1.5 False Alarm Test. — Structural response can change due to operational/environmental variability 
or increase in stiffness. To study the false alarm, two locations on a simply supported plate are 
considered. At first, both locations are weakened by decreasing the stiffness and then strengthened 
another time by increasing the stiffness. Damage energy vectors are extracted to distinguish the false 
alarm case with the actual damage case. Here, the damage vectors will point in the direction of decreased 
total energy. The vectors pointing outward from a location indicates weakening of the material, whereas, 
the vectors pointing inward at a location indicates strengthening of the material. 

5.1.6 A comparative case study. — J.-C. Hong, et al. (ref. 48) described that wavelet transforms can be 
adopted for damage detection using the vibration modes of a beam. They also used the wavelet toolbox in 
Matlab for the simulations. A comparative case study is performed with the simulations described in 
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reference 48 using the wavelet toolbox in Matlab. In this section of the study, both fine and coarse 
networks are used to observe the behavior of wavelet analysis in global damage detection. 

The fine network is used here for comparison purpose only and also because it is the only case 
presented in the above reference 48. However, from a practical viewpoint, coarse networks are crucial in 
judging the effectiveness of any scheme. These coarse network simulations are performed here only using 
the detection scheme proposed in this report, since unfortunately there were no comparable studies in the 
available literature searched using other global detection schemes. 

Along with these simulation and comparative studies, some actual experiments, which are provided 
by NASA Glenn Research Center, were also analyzed. These experiments are conducted by Electronic 
Speckle Pattern Interferometry (ESPI) technique, which is used to determine the correct overall 
underlying mode shapes and frequencies. Actually, these experimental analyses motivated this study 
towards investigating the effect of noisy measurements and the feature reconstruction. The results of 
above cases are presented in the following sections. 


5.2 Simulation Results 


5.2.1 Low Cyclic Fatigue Failure Mode . — Damage in a structure occurs after a low number of cyclic 
loadings in Low Cyclic Fatigue failure mode. This type of damage is simulated by a sudden degradation 
of the Modulus of Elasticity of element/elements in the finite element model. A metallic plate shown in 
figure 5.1 is considered for the simulations with the damage oriented horizontally. 


7 in. 




5.25 in. 


Material properties: 

Mass of the plate, m: 

Poison’s ratio, v : 

Modulus of elasticity, E\ 

Modulus of elasticity of damaged elements, E d \ 
Percentage of material degradation: 

Percentage of damaged area: 

Location of the damage from left top comer: 


7.67E-04 lb 
0.3 

2.96E+07 ksi 

1.00936E+07 ksi 

66 percent of the actual Modulus 

0.17 percent of total area 

2.0625 in. in v and 1.5 in. in y — directions 
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A metallic plate with the above material properties are considered for the simulations. The plate is 
fixed at the bottom edge and left free at other edges. A mesh of size 40 by 40 is used. A total of six 
elements are damaged, i.e., the Modulus of Elasticity is reduced to E d from E , at a location treated as a 
single damage. Figure 5.2(a) shows the first five mode shapes of the cantilever plate and they can be 
described as first bending mode, first twisting mode, first mixed (bending and twist) mode, second 
bending mode, and second twisting mode, respectively. Figure 5.2(b) shows the damage parameter 
contours in which we can view the location and size of the damage. Figure 5.3 shows the shapes of 
undamaged, damaged and delta (undamaged - damaged) mode shapes. We can see that the mode shapes 
of damaged and undamaged are similar, but delta mode shapes are completely different than the actual 
mode shape of the plate. Figure 5.4 shows the damage energy vectors. The vectors point in the direction 
of decreased total energy. In all the mode shapes we can see that the vectors are going away from a 
particular location. That means, that location is weaker compared to the other portions of the plate. 

The undamaged and the damaged frequencies are tabulated in tables 5.1 and 5.2 respectively. When 
we compare them, there is less than 0.01 percent reduction in the damaged frequencies in first five modes. 
Here, the frequencies do not differ much with the damage. In other words, we can say that, the damage is 
very tiny since it does not significantly change the frequencies in the damaged plate. This is also further 
confirmed if one compare the very small difference in the mode changes of figure 5.3(c), as compared to 
the rather large magnitudes in the almost-identical overall underlying modes (before and after damages) 
in figures 5.3(a) and 5.3(b). This is of course of practical significance in that the detection at early stages 
of small-size defects is always desirable. Furthermore, experimental methods needed for the detection are 
distinctly more demanding in their level of accuracy compared to those that are merely targeting the 
overall mode shapes as in the traditional modal testing. 


TABLE 5.1.— UNDAMAGED FREQUENCIES (CYC/SEC) OF 
MATERIAL DEGRADATION - SINGLE DAMAGE 


Mode-1 

149.2480086 

Mode-2 

295.1841725 

Mode-3 

722.7485231 

Mode-4 

942.1421981 

Mode-5 

1157.100669 


TABLE 5.2.— DAMAGED FREQUENCIES (CYC/SEC) OF 
MATERIAL DEGRADATION - SINGLE DAMAGE 


Mode-1 

149.2214258 

Mode-2 

294.7049912 

Mode-3 

721.222923 

Mode-4 

939.1733984 

Mode-5 

1155.178122 
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1 0.417 
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Mode -5 



39.5 

36.2 
32.9 

29.6 

26.3 

23.1 
19.8 
16.5 

13.2 
9.90 
6.61 
3.31 

0.0244 


Figure 


(a) (b) 

5.2. — Metallic Plate with Single Damage (a) Mode Shapes (b) 


Damage Parameter 
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Mode-1 
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Mode-2 
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Mode-4 
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0.0163 
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I 0.00786 
I 0.00576 
I 0.00366 
I 0.00156 
1-0.000535 
| -0.00263 
I -0.00473 
* -0.00683 



Mode-5 



1 0.00602 
0.00491 
0.00380 
0.00269 
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0.000474 
-0.000636 
•0.00175 
•0.00286 
-0.00397 
-0.00508 
-0.00619 
■0.00730 


(a) (b) (c) 

Figure 5.3. — Metallic plate with single damage (a) Undamaged mode shape (b) Damaged 
mode shape (c) Delta mode shape (undamaged shape-damaged mode shape). 
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Figure 5.4. — Metallic Plate with Single Damage - Damage nodal force vectors 
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////////////////// 


Fixed Edge 

Figure 5.5. — Cantilever plate with multiple damages. 


In the following case, two damages are simulated as in figure 5.5. The material properties and the 
dimensions for damage 1 have taken from the previous case. 

For the damage2 : 

Modulus of elasticity of damaged elements for damage2, E d2 : 1.04488E+07 ksi 

Percentage of material degradation: 65 percent of the actual Modulus 

Percentage of damaged area: 0.40 percent of total area 

Location of the damage2 from left top comer: 5.075 in. in x and 3.625 in. in y-directions 

Here, seven elements are damaged, i.e., Modulus of Elasticity is reduced to E d2 from E, at the location 
of damage2. Figures 5.6(b) and 5.8 show the damage parameter contours and damage vectors of first five 
modes respectively. The damage can be identified by observing the first five modes. Here, we can not 
locate both damages in all five modes, but we can detect the damage if one can provide the first few 
fundamental modes. Therefore, we can conclude that we need at least a few modes to confirm the 
existence and location of the damage. Also we can observe that delta mode shapes in figure 5.7(c) are 
completely different compared to the single damage case. When we compare the undamaged and the 
damaged frequencies in tables 5.3 and 5.4, damaged frequencies are reduced by less than 0.01 percent in 
first five modes. The same behavior was observed in single damage case also. As in the single damage 
case, the difference in the mode changes is very small. Therefore, we can conclude that a tiny amount of 
damages in structure will not change much in the frequencies and mode shapes, and the level of accuracy 
needed in experimental methods in the detection is more demanding issue. 


TABLE 5.3.— UNDAMAGED FREQUENCIES (CYC/SEC) OF 
MATERIAL DEGRADATION - MULTIPLE DAMAGES 


Mode-1 

149.2480086 

Mode-2 

295.1841725 

Mode-3 

722.7485231 

Mode-4 

942.1421981 

Mode-5 

1157.100669 


TABLE 5.4.— DAMAGED FREQUENCIES (CYC/SEC) 
OF MATERIAL DEGRADATION - MULTIPLE DAMAGES 


Mode-1 

148.9227814 

Mode-2 

294.2309847 

Mode-3 

720.2086541 

Mode-4 

937.8735199 

Mode-5 

1154.13901 
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(b) 


Figure 5.6. — Metallic plate with two damages (a) Mode shapes (b) Damage parameter. 
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(a)(b) (c) 

Figure 5.7. — Metallic plate with two damages (a) Undamaged mode shape (b) Damaged mode shape 
(c) Delta mode shape (undamaged shape-damaged mode shape). 
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Figure 5.8. — Metallic plate with two damages - damage nodal force vectors. 
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5.2.2 High Cyclic Fatigue Failure Mode . — The damage in a structure can occurs after a high number 
of cyclic loadings in high Cyclic Fatigue failure mode. This type of damage/line crack is simulated by 
creating multiple nodes at the same geometric location in the finite element model. A metallic plate in 
figure 5.9 is considered for simulations. The plate material properties are taken same as previous cases. 

A crack of length 175 mil (1 mil = 0.001 in.), oriented horizontally, on a plate fixed at one edge is 
simulated with a mesh of 20 by 10 size. Initially a mesh of 80 by 20 size is generated and an extra node is 
created at the location of damage. Nodal connectivity is changed using the extra node to simulate the 
crack. Then the equivalent nodal displacements on a mesh of 20 by 10 size are considered to extract the 
first five fundamental modes. The results are presented from figures 5.10 to 5.12. In figure 5.12, damage 
vectors of mode 5 are shown at two locations unlike other modes. This is happened due to the false alarm, 
which is caused by the processing of huge data files. But, overall we can confirm the existence and 
location of the damage by considering the first five fundamental modes. Therefore, to confirm the 
existence and the location of the damage we require at least few mode shapes. The table 5.5 shows the 
undamaged frequencies of the first five modes. 


TABLE 5.5.— UNDAMAGED FREQUENCIES (CYC/SEC) OF LINE CRACK 


Mode-1 

149.1728511 

Mode-2 

294.822034 

Mode-3 

721.1457221 

Mode-4 

943.1049457 

Mode-5 

1158.669922 
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Figure 5.9. — Cantilever plate with line crack. 
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(a) (b) 

Figure 5.10. — Metallic plate with a line crack (a) Mode shapes (b) Damage parameter. 
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Figure 5.11. — Metallic plate with a line crack (a) Undamaged mode shape 
(b) Damaged mode shape (c) Delta mode shape. 
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Figure 5.12. — Metallic plate with a line crack - damage nodal force vectors. 
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5.2.3 Corrosion in the Structure . — Corrosion of the material is also a type of damage to the 
structure. When corrosion of the structure begins, we can observe tiny holes. This is simulated by eroding 
an element from the mesh geometry. Two different mesh sizes, 40 by 40 and 120 by 120, are considered 
for this case. The plate in figure 5.13 is considered with the material properties as in previous cases. 

A finite element model is created using ABAQUS with different mesh sizes and an element is deleted 
from the geometry to simulate the corrosion effect. The first five fundamental modes are generated using 
frequency analysis and the undamaged frequencies are listed in the tables 5.6 and 5.7. Displacements of 
each mode are saved and entered into the data input file of the detection scheme for damage detection. 

The features presented in the figures from figures 5.14 to 5.19 describe the existence and the location 
of the damage in both cases. The figures 5.14 and 5.16 show the existence of the damage for 40 by 40 
mesh analysis, where 0.0625 percent of total area is damaged and the figures 5.17 and 5.19 show the 
existence of the damage for 120 by 120 mesh analysis, where 0.0067 percent of total area is damaged. As 
expected, in both cases, the delta mode shapes, presented in figures 5.15 and 5.18, are completely 
different than the actual mode shapes. 


TABLE 5.6.— UNDAMAGED FREQUENCIES (CYC/SEC) 
OF CORROSION - MESH 40 BY 40 


Mode-1 

144.133321 

Mode-2 

285.0233531 

Mode-3 

698.1044529 

Mode-4 

910.8488979 

Mode-5 

1118.65208 


TABLE 5.7.— UNDAMAGED FREQUENCIES (CYC/SEC) 


OF CORROSION - MESH 120 BY 120 


Mode-1 

149.2513612 

Mode-2 

295.1571222 

Mode-3 

722.5704527 

Mode-4 

941.9485506 

Mode-5 

1157.162608 
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k ►! 



Location of the damage from left top comer: 

Percentage of damaged portion (40 by 40 mesh): 
Percentage of damaged portion (120 by 120 mesh): 


2.0625 in. in x - direction 
1 .5 in. in y - direction 
0.0625 percent of total area 
0.0067 percent of total area 
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Figure 5.14. — Metallic plate (40 by 40 mesh) with a tiny hole due 
to corrosion (a) Mode shapes (b) Damage parameter. 
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Figure 5.15. — Metallic plate (40 by 40 mesh) with a tiny hole due to corrosion 
(a) Undamaged mode shape (b) Damaged mode shape (c) Delta mode shape. 
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Figure 5.16. — Metallic plate (40 by 40mesh) with a tiny hole-damage nodal force vectors 
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(a) (b) 

Figure 5.17. — Metallic plate (120 by 120 mesh) with a tiny hole 
due to corrosion (a) Mode shapes (b) Damage parameter. 
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Figure 5.18. — Metallic plate (120 by 120 mesh) with a tiny hole due to corrosion 
(a) Undamaged mode shape (b) Damaged mode shape (c) Delta mode shape. 
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Figure 5.19. — Metallic plate (120 by 120 mesh) with a tiny hole-damage nodal force vectors. 


5.2.4 Sensor Network on the Structure in Global Damage Detection . — The effect of the sensor 
network is a major issue in global damage detection techniques. Therefore, a case study is performed on 
the density of sensor network. Three sizes of arrangement of sensor networks are considered. For 
simulations, three mesh sizes are taken viz. extensive network (40 by 40), intermediate network (20 by 
20), and coarse network (10 by 10). The case in section 5.2.1 is considered as an extensive network and 
equivalent displacements of coarser networks are extracted from this extensive network. 
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Important features are extracted from the simulations using Eview. Mode shapes, delta mode shapes, 
and damage parameter values are extracted for comparison. In all three cases, mode shapes are the same 
as in cantilever plate. As expected, delta mode shapes are completely different than the actual mode 
shapes. We are able to extract the feature of the damage for every size of the network arrangement. If we 
observe the results carefully, the size of the damage is widened from fine to coarse network. Since this is 
a finite element analysis, extraction of size of the damage is proportional to the size of the elements in the 
mesh. 

The proposed detection scheme in this report is robust in the extraction of the damage even with 
coarser sensor networks. The required figures are presented from the figures 5.20 to 5.22. 
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Figure 5.20. — Mode shapes (a) Extensive sensor network 
(b) Intermediate sensor network (c) Coarse sensor network. 
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Figure 5.21. — Delta mode shapes (a) Extensive sensor network 
(b) Intermediate sensor network (c) Coarse sensor network. 
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Figure 5.22. — Damage parameter (a) Extensive sensor network 
(b) Intermediate sensor network (c) Coarse sensor network. 

5 . 2.5 Analysis of Experimental Files . — -NASA Glenn Research Center is using Electronic Speckle 
Pattern Interferometry (ESPI) experimental setup for the measurement of overall underlying mode shapes, 
which is used here for the Global Damage Detection. Experimental files provided by NASA are analyzed 
to detect the presence of the damage. When we observe the undamaged, damaged, and delta mode shapes, 
the delta mode shapes have similar pattern to that of the actual mode shapes which is unlike to the 
previous cases. The damage parameter also shows the same pattern, instead of showing sharp damage. 

The undamaged frequencies are presented in table 5.8. The figures 5.23 and 5.24 show the mode 
shapes and damage parameter contours respectively. 
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TABLE 5.8.— UNDAMAGED FREQUENCIES (CYC/SEC) 
OF EXPERIMENTAL RESULTS. 
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Figure 5.23. — Experimental results with single damage (a) Undamaged mode shape 
(b) Damaged mode shape (c) Delta mode shape. 
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Figure 5.24. — Experimental results with single damage - damage parameter. 


5.2.6 Effect of the Noise in Measurements .—To study the effect of the noise in the detection 
scheme, white noise and pink noise are added to first five fundamental mode shapes. The damage 
parameter is shown to be able to extract the damage. In addition, each mode is able to withstand a certain 
level of noise. The ideal (noise free) case is considered from section 5.2.1. 

From figures 5.25 to 5.29 shows the damage parameter contours for fist five modes. In particular, 
figure 5.25(b) shows the damage parameter of mode-1 with white noise as vectors, where we can increase 
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the vector threshold to retain those vectors with high magnitudes. Each figure shows the damage 
parameter contours of a mode with various levels of noise, which include both white and pink noises. 
These figures emphasize that each mode have different critical levels of noise. The level of noise that 
each mode can withstand is summarized in table 5.1. 
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Figure 5.25. — Mode-1 with white noise (a) damage parameter contour view. 
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Figure 5.25. — Mode-1 with white noise (b) damage parameter vectors view. 
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Figure 5.26. — Mode-1 with pink noise 
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Figure 5.27. — Mode-2 with white noise 
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Figure 5.28. — Mode-2 with pink noise 
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Figure 5.29. — Mode-3 with white noise and pink noise. 
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Figure 5.30. — Mode-4 with white noise and pink noise. 
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Figure 5.31. — Mode-5 with white noise and pink noise. 
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TABLE 5.9.— ROBUSTNESS OF GLOBAL DAMAGE DETECTION WITH WHITE NOISE AND PINK NOISE 


Out-of-plane 

Mode-1 

Mode-2 

Mode-3 

Mode-4 

Mode-5 

White noise level 

Up to 175% 

Up to 105% 

Up to 17.5% 

Up to 70% 

Up to 70% 

Pink noise level 

Up to 175% 

Up to 175% 

Up to 35% 

Up to 105% 

Up to 105% 


From table 5.9, we can say that, if one can provide several modes, this detection scheme can able to 
withstand up to 70 percent of noise. 

5.2. 7 Signal - Feature Reconstruction . — To study signal feature reconstruction, white noise is added 
to the measurements in two ways. 

1. ) Noise is added to the displacements in order to see its effect on the damage parameter and the 
subsequent reconstruction of the damage parameter feature using Wavelet Toolbox in Matlab. 

2. ) Noise is added to the damage parameter values and the reconstruction damage parameter feature 
using Wavelet Toolbox in Matlab. The following equation is used to add the noise in the damage 
parameter values. 


P n = P + max(P)*% of N 


where, 

P n = Noised damaged Parameter 
P = Ideal Damage Paramter 
N = Noise = rand(-l, 1) 

Ideal (noise free) case is considered from the section 5.2.1. Results of case-1, which is the addition of 
the noise to the displacements, simulations are presented in figure 5.32 and results of case-2, which is the 
addition of the noise to the damage parameter values, simulations are presented in figure 5.33. 
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Figure 5.32.— Reconstruction of damage parameter 
with white noise added to the displacements. 
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Figure 5.33. — Reconstruction of damage parameter with 
white noise added to the damage parameter values. 
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An additional case is considered which is unlike to the previous cases. Here, instead of taking a single 
damage, two damages are considered on a plate shown in figure 5.34, one is oriented horizontally and 
other is oriented at 45°. The plate is simply supported and it has the same dimensions and material 
properties as in the section 5.2.1. Modulus of elasticity is reduced by 66 percent for both of the damages. 
The five fundamental mode shapes are extracted from the dynamic analysis and the corresponding 
undamaged frequencies are listed in table 5.10. The white noise is added to the damage parameter to 
study the feature reconstruction capability in global damage detection. Mode shapes and damage 
parameter contours are presented in figure 5.35. Figure 5.36 shows the noised and the reconstructed 
damage parameter contours. 
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h ►! 



0.655 in. \ 
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\ 0.525 in. 
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0. 22 in. | 


::: 0.0625 in. damage2 


Figure 5.34 — Simply supported plate with two damages. 


Location of damagel from left bottom corner: 1.05 in. in v-dir and 1.05 in. iny-dir 
Location of damage2 from right top corner: 1.225 in. in v-dir and 1.31 in. iny-dir 
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Figure 5.35. — Simply supported plate with two damages oriented horizontally and inclined at 45° 
(a) Mode shapes (b) Damage parameter (noise free). 
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Figure 5.36. — Reconstruction of damage parameter with white noise added to damage 
parameter values with two damages oriented horizontally and inclined at 45°. 
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TABLE 5.10.— UNDAMAGED FREQUENCIES (CYC/SEC) 
OF SIMPLY SUPPORTED PLATE 


Mode-1 

231.0301748 

Mode-2 

473.1115781 

Mode-3 

934.7355844 

Mode-4 

1245.773229 

Mode-5 

1316.060446 


5.2.8 False Alarm Test . — A false alarm case study is performed to investigate the proposed detection 
scheme for this research. A simply supported plate of 7 in. by 5.25 in. by 0.125 in./is considered. Two 
areas are considered for damage, one is oriented horizontally and other is oriented at 45°. Figure 5.34 
shows the size and locations of the damages. For the false alarm test, at first the structure is weakened by 
66 percent and in the second test the structure is strengthened by 66 percent. Damage vector fields are 
extracted to distinguish the weakening and strengthening of the structure at a specific location. The 
damage vector points in the direction of decreased total energy. Figure 5.52 shows the weakening of the 
structure, where the damage vectors are directed outward from the location. Figure 5.53 shows the 
strengthening of the structure, where the damage vectors are directed inward at the location. 
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Figure 5.37. — Damage nodal force vectors for weakening of the structure. 
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Figure 5.38. — Damage nodal force vectors for strengthening of the structure. 
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5.2.9 A Comparative Case Study. — J.-C. Hong, et al. (ref. 48) described how wavelet transforms can 
be adopted for damage detection using the vibration modes of beam. They suggested Continuous Wavelet 
Transforms (CWT) with the Mexican-hat wavelet found in the Matlab wavelet toolbox. The Mexican-hat 
wavelet function is the second derivative of the Gaussian probability density function. They also used 
ANSYS to simulate a simply supported beam model with 2400 two-node beam elements (element size = 
0.5 mm). Location of the damage is at a = 800 mm. Material properties are Young’s Modulus = 70 GPa, 
Density = 2700 kg/m 3 . 

The model for the analysis is shown in figure 5.37. Modal analysis is performed in ANSYS and the 
first bending mode is passed through the wavelet toolbox. The wavelet analysis has detected the damage 
and has located it exactly at the correct position. Figure 5.39 shows the contour plot of the wavelet 
transform in which we can observe the location of the damage. 

A similar case study is performed with the simulations described in the literature using ABAQUS. figure 
5.41, which is a contour of wavelet analysis, shows the location of the damage. Here, the wavelet analysis 
utilized an extensive (2400 points) sensor network. As we can clearly observe, there is a measurement 
(sensor) right on top of the damage. In reality, there is a small chance that there will be a measurement 
right on top of the damage since location of the damage is unknown and arrangement of extensive sensor 
networks also very difficult. Therefore, coarser network is considered. Instead of passing 2400 data 
points, in first case, only 100 points are passed (a measurement after every 24 points from 2400 data 
points) and in second case 160 data pointes are passed (a measurement after every 15 points from 
2400 data points). In both cases, the direct wavelet approach is not able to detect the damage. The results 
are presented form figures 5.42 to 5.45. 
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w 




_! 

Isj 


( L=120D mm, h=20min. w =0.5 mm) 
Figure 5.39. — Model for the simulation. 


NAS A/CR— 2006-2 1 4225 


93 




Amplitude 



f 'rad; luatloi 


l 


/ H 


Jf \ 


/ A 

. 

/ \ 

/ 

\ 

\ 

/ 

/ 

Btflcii vr 1 iliciui a c rack 1 

/ 

Riis*n ih Nil si crt«' k 

/ 

\ 

■ / 

\ 

\ 

/ 


/ 

\ 

, / 

\ . 

/ 

\ 


0 240 440 604 1*00 1444 1244 


fa) Lofalioi i (mm> 

Figure 5.40. — First bending 
mode from the literature. 



m 400 m m logo uoo 

(b) Location (mn) 


Figure 5.41 . — Contour plot of wavelet 
transform from the literature. 



Figure 5.42. — First bending mode. 



Figure 5.43. — Contour plot of wavelet transform. 
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Figure 5.45. — Contour plot of wavelet transform with 100 data points. 



Figure 5.46. — First bending mode with 160 data points. 



Figure 5.47. — Contour plot of wavelet transform with 160 data points. 


An equivalent simulation is performed using the damage detection scheme proposed in this report for 
a simply supported beam of same dimensions. Damage is simulated over a 0.1 percent area of the beam 
by degradation of the modulus of elasticity by 14 percent at the location of the damage. Instead of passing 
the mode shape to the wavelet analysis, the damage parameter is passed which is derived from the 
proposed scheme. Figure 5.47 shows the contour of wavelet analysis. Length / is the approximate size of 
the damage shown in the contour with the displacements measured at 600 points. In addition, the damage 
parameter is calculated for coarser networks also. The displacements are considered at 20 and 10 points to 
simulate the coarser networks. Care was taken such that there is no measurement at the location of the 
damage. In both coarser networks, the damage was successfully detected. If we observe the wavelet 
analysis contour plots in figures 5.49 and 5.51, the size of the damage is wider than in the figure 5.47, 
since it depends on the element size. 
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Figure 5.48 Damage Parameter in first bending mode with 600 data points. 
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Figure 5.49. — Contour plot of wavelet transform with 600 data points. 



Figure 5.50.— Damage Parameter in first bending mode with 20 data points. 



Figure 5.51. — Contour plot of wavelet transform with 20 data points. 
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Figure 5.52. — Damage Parameter in first bending mode with 10 data points. 



◄ ► 

12 

Figure 5.53. — Contour plot of wavelet transform with 10 data points. 
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Figure 5.54. — Cantilever plate with 12 damaged elements. 
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5.2.10 A Case Study on the Severity of the Damage . — A robust damage detection technique should 
be able to judge the severity of the damage. To study the severity of the damage two different damage 
sizes are considered. One size of the damage has considered from the section 5.2.1. To simulate more 
severe damage, double the size of the damage in the section 5.2.1 is considered. Including the six damage 
elements in section 5.2.1, six more elements are damaged. Figure 5.52 shows the size and shape of the 
damage. The mode 3, which is the mixture of bending and twisting in cantilever plate, is considered to 
investigate the severity. 

Figure 5.53 shows the delta mode shape, which the difference between the damaged and the 
undamaged mode shapes, wherein the difference is greater in the double the size of the damage than the 
actual size. Figure 5.54 shows the damage parameter contour, wherein the magnitude of the double the 
size of the damage is greater then the actual size of the damage. The two figures show the severity of the 
damage. Therefore, we can judge the severity of the damage by inspecting the delta mode shape and the 
damage parameter contours. 
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Figure 5.55. — Delta mode shapes (a) Damage with actual size (b) Damage with double the size. 
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Figure 5.56. — Damage parameter (a) Damage with actual size (b) Damage with double the size. 
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Chapter VI 

Summary, Conclusions and Future Work 

6.1 Summary 

The presented study was conducted in conjunction with an ongoing research program dealing with the 
development of an integrated analytical/experimental nondestructive evaluation (NDE) methodology for 
global damage detection. The important ingredient is the use of a damage energy parameter. It provides 
the required sensitivity measure for detection of the damage on the basis of measured structural response. 

From the viewpoint of practical utilization in structures, extensive testing for a valid assessment of 
any NDE technique becomes necessary. Especially in view of the many complicating factors and the 
variety of damage scenarios that are likely to be encountered, assessment of the scheme becomes a 
necessary requirement. In the real time environment, damage can manifest itself as: weakened material, a 
line crack, or corrosion. An investigation is performed keeping the above types of damage in mind. A 
case study is conducted on the density of the arrangement of the sensor network. The extensive/fine 
sensor network has unlimited capability of measuring response at all regions in the structure. However, 
from a practical viewpoint, it is very difficult to arrange the extensive network due to various factors. 
Therefore, effectiveness of the scheme is tested with coarse sensor networks. In addition, by observing the 
difficulty in obtaining noise free measurements in the real time environment, a case study is performed to 
assess the capability of the scheme to detect damage at various noise levels. If the technique is 
excessively sensitive to the noise levels then the effect of the “true” damage will be completely 
“overshadowed” by the noise. That means, extracting the damage feature will be a difficult task. 
Therefore, feature reconstruction plays an important role in global damage detection. In this report, the 
feature extraction capability of the detection scheme is investigated extensively. With an eye towards 
reducing the noise intensity, the development of different types of filters has been undertaken. Here, 
wavelet theory is utilized to reconstruct the damage feature from the noise to enhance the feature 
extraction. An important factor in global damage detection is to over come the false alarm, i.e., the 
detection scheme should able to detect whether there is a weakening or a strengthening of the structure. 

As a result, a false alarm test is conducted on the detection scheme. 


6.2 Conclusions 

Based on the presented results, the following conclusions can be made: 

1. The detection scheme is able to detect and locate any type of damage (i.e., the different failure 
modes) present in the structure. As expected, the so-called delta mode shapes are completely different 
than the actual mode shapes. 

2. The required features are extracted efficiently by the Eview graphical program, which allows the 
user to view the large and sometimes complex files generated by the proposed detection scheme. 

3. The detection scheme is efficient enough to extract the feature of the damage even with coarse 
sensor networks. 

4. When actual preliminary experimental files (using Electronic Speckle Pattern Interferometry 
(ESPI)) were analyzed, the delta mode shapes behaved similar to the actual mode shapes. These 
experimental measurements were determined not to be directly usable for the purpose of detecting small 
amounts of damages. As they presently exist, the ESPI techniques focus on the determination of correct 
overall underlying mode shapes and frequencies. Recall that for the small amounts of damages targeted in 
the present study, essentially identical overall modes of vibration were obtained before and after damage 
(see figs. 5.3(a) and (b)).This indeed is one of the main motivations for conducting the noise level 
quantification studies reported here. 
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5. When the noise is added to the displacements, i.e., before processing the data, the damage 
detection scheme has been shown to be able to withstand up to 70 percent of noise (see table 5.9) when 
several modes are provided for detection. 

6. When noise is added to the damage parameter, i.e., after processing the data, the damage 
detection scheme has shown to be able to withstand up to 50 percent of noise (see figs. 5.33 and 5.36). 

7. A one dimensional wavelet analysis is used for the feature reconstruction of the damage 
parameter signal from the noise. Wavelet transforms greatly reduced the intensity of the noise by 
retaining the significant damage values. 

8. The detection scheme is able to distinguish the strengthening and the weakening of the structure, 

i.e., it overcame the false alarm test. 

9. We can conclude that, if ideal measurements are provided, even though there may be few 
measurements, i.e., on coarser sensor network, the proposed scheme can detect and locate the damage. 
Also, the scheme can withstand up to a level of 70 percent of noise when raw displacements are directly 
utilized. 

10. However, utilizing further enhancement through wavelet tools, i.e., denoising/ feature 
reconstruction, one is able to use the detection scheme with 105 percent of noise in the displacements 
when several modes are provided. 

1 1 . Enhancement of the damage feature extraction capability by reducing the intensity of the noise is 
achieved using the wavelets (see figs. 5.32, 5.33, and 5.36). 
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Appendix A 

C++ Programs for Global Damage Detection 


A.l Program to Add the Noise to the Displacements in the Data Input File 


///////////PROGRAM TO ADD THE NOISE TO THE DATA INPUT (.din) FILE ///////// 

#include<iostream> 

#include<iomanip> 

#include<cstdlib> 

#include<fstream> 

#include<math.h> 

using namespace std; 

//int main() 

int main(int argc ,char * argv[]) 

{ 

if(argc != 5) 

{ 

cout«"Error 

cout«"Usage: commandline< inputfile noise 1 file noise 1 file result file >"«endl; 

return - 1 ; 

} 

char buffer[2093]; 
int intTmp; 
int maxnode; 
int maxelement; 
double m,ml,m2; 

// int count=0; 

//ifstream inFile_din( M phasel_format.din M ,ios::in); 
ifstream inFile_din(argv[ 1 ],ios: :in); 

if(!inFile_din) 

{ 

cerr«". din file not found! !!!!!! M «endl; 

} 

inFile_din.getline(buffer,2093); 

inF ile_din»intT mp ; 

inF ile_din»max_node; 

inFile_din.getline(buffer,2093); 

inFile_din.getline(buffer,2093); 

inFile_din.getline(buffer,2093); 
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while(l) 

{ 

inF ile_din»intT mp ; 
inFile_din.getline(buffer,2093); 
inFile_din.getline(buffer,2093); 
if(intTmp == max_node) 
break; 

} 

while(l) 

{ 

inF ile_din»intT mp ; 
inFile_din.getline(buffer,2093); 


} 


if(intTmp == maxnode) 
break; 


inFile_din.getline(buffer,2093); 

inF ile_din»intT mp ; 

inF ile_din»max_element ; 


int mat; 

inF ile_din»mat ; 
inFile_din.getline(buffer,2093); 

for(int i=0; i<mat; i++) 

{ 

inFile_din.getline(buffer,2093); 

} 

while(l) 

{ 

inF ile_din»intT mp ; 

inFile_din.getline(buffer,2093); 

inFile_din.getline(buffer,2093); 

if(intTmp == maxelement) 
break; 

} 

inFile_din.getline(buffer,2093); 


// Calculating Norm of UnDamaged Displacements 

float Unorml, Unorm2, Unorm3, Unorm4, Unorm5, maxUnorm = 0.0; 
Unorml = Unorm2 = Unorm3 = Unorm4 = Unorm5 = 0.0; 
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bool done = false; 
int x=0; 

double U1 [10000], U2[10000],U3[10000],U4[10000],U5[10000]; 

while(ldone) 

{ 

inF ile_din»intT mp ; 
if(intTmp == max_node) 
done = true; 

inF ile_din»intT mp ; 
inF ile_din»intT mp ; 

inFile_din»ml ; 

U 1 [x]=ml ; //storing the undamaged x-displacements 
Unorml = Unorml+ml*ml; 

inF ile_din»intT mp ; 
inFile_din»ml; 

U2[x]=ml ; //storing the undamaged y-displacements 
Unorm2 = Unorm2+ml*ml; 

inF ile_din»intT mp ; 
inFile_din»ml ; 

U3 [x]=ml ; //storing the undamaged z-displacements 
Unorm3 = Unorm3+ml*ml; 

inF ile_din»intT mp ; 
inFile_din»ml; 

U4[x]=ml ; //storing the undamaged Rx-displacements 
Unorm4 = Unorm4+ml *ml ; 

inF ile_din»intT mp ; 
inFile_din»ml ; 

U5 [x]=ml ; //storing the undamaged Ry-displacements 
Unorm5 = Unorm5+ml*ml; 

x++; 

} 

double udSignall = Unorml; 
double udSignal2 = Unorm2; 
double udSignaB = Unorm3; 
double udSignal4 = Unorm4; 
double udSignal5 = Unorm5; 
double udSignal; 

udSignal = Unorml +Unorm2+Unorm3+Unorm4+Unorm5; 
Unorml = sqrt(Unorml); 
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Unorm2 = sqrt(Unorm2); 

Unorm3 = sqrt(Unorm3); 

Unorm4 = sqrt(Unorm4); 

Unorm5 = sqrt(Unorm5); 

inFile_din.getline(buffer,2093); 

inFile_din.getline(buffer,2093); 

inFile_din.getline(buffer,2093); 

// Calculationg Norm of Damaged Displacements and Norm of Delta 

float DeltaNorml ,Dnorml ,Dratio 1 ,Uratio 1 ,DeltaNorm2,Dnorm2,Dratio2,Uratio2, 
DeltaNorm3 ,Dnorm3 ,Dratio3 ,Uratio3 , DeltaNorm4,Dnorm4,Dratio4,Uratio4, 
DeltaNorm5,Dnorm5,Dratio5,Uratio5, maxDnorm = 0.0; 

double D1 [10000], D2[10000],D3[10000],D4[10000],D5[10000]; 
double Delta; 

done = false; 
x=0; 

DeltaNorml =Dnorml = 0.0; 

DeltaNorm2=Dnorm2 = 0.0; 

DeltaNorm3=Dnorm3 = 0.0; 

DeltaNorm4=Dnorm4 = 0.0; 

DeltaNorm5=Dnorm5 = 0.0; 
while(ldone) 

{ 

inF ile_din»intT mp ; 
if(intTmp == max_node) 
done = true; 

inF ile_din»intT mp ; 
inF ile_din»intT mp ; 

inFile_din»ml ; 

Dnorml = Dnorml+ ml*ml; 

Dl[x] = ml; 

Delta=D 1 [x] -U 1 [x] ; //(Damage - undamage = Delta) 

DeltaNorml = DeltaNorml+Delta*Delta; 

inF ile_din»intT mp ; 
inFile_din»ml; 

Dnorm2 = Dnorm2+ ml*ml; 

D2[x] = ml; 

Delta=D2[x]-U2[x]; //(Damage - undamage = Delta) 

DeltaNorm2 = DeltaNorm2+Delta*Delta; 

inF ile_din»intT mp ; 
inFile_din»ml; 
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Dnorm3 = Dnorm3+ ml*ml; 

D3[x] = ml; 

Delta=D3[x]-U3[x]; //(Damage - undamage = Delta) 
DeltaNorm3 = DeltaNorm3+Delta*Delta; 

inF ile_din»intT mp ; 
inFile_din»ml ; 

Dnorm4 = Dnorm4+ ml*ml; 

D4[x] = ml; 

Delta=D4[x] - U4[x]; //(Damage - undamage = Delta) 
DeltaNorm4 = DeltaNorm4+Delta*Delta; 

inF ile_din»intT mp ; 
inFile_din»ml ; 

Dnorm5 = Dnorm5+ ml*ml; 

D5[x] = ml; 

Delta=D5[x]- U5[x]; //(Damage - undamage = Delta) 
DeltaNorm5 = DeltaNorm5+Delta*Delta; 

x++; 

} 

double dSignall = Dnorml; 
double dSignal2 = Dnorm2; 
double dSignaB = Dnorm3; 
double dSignal4 = Dnorm4; 
double dSignal5 = Dnorm5; 
double dSignal; 

dSignal = Dnorml +Dnorm2+Dnorm3+Dnorm4+Dnorm5; 

Dnorml = sqrt(Dnorml); 

Dnorm2 = sqrt(Dnorm2); 

Dnorm3 = sqrt(Dnorm3); 

Dnorm4 = sqrt(Dnorm4); 

Dnorm5 = sqrt(Dnorm5); 

DeltaNorml = sqrt(DeltaNorml); 

DeltaNorm2 = sqrt(DeltaNorm2); 

DeltaNorm3 = sqrt(DeltaNorm3); 

DeltaNorm4 = sqrt(DeltaNorm4); 

DeltaNorm5 = sqrt(DeltaNorm5); 

if(Unorml=0) Uratiol =0; 

else Uratiol = DeltaNorml /Unorml; 

if(Unorm2==0) Uratio2 =0; 
else Uratio2 = DeltaNorm2/Unorm2; 

if(Unorm3==0) Uratio3 =0; 

else Uratio3 = DeltaNorm3/Unorm3; 
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if(Unorm4==0) Uratio4 =0; 

else Uratio4 = DeltaNorm4/Unorm4; 

if(Unorm5==0) Uratio5 =0; 

else Uratio5 = DeltaNorm5/Unorm5; 

if(Dnorml==0) Dratiol =0; 

else Dratiol = DeltaNorml/Dnorml; 

if(Dnorm2==0) Dratio2 =0; 

else Dratio2 = DeltaNorm2/Dnorm2; 

if(Dnorm3=0) Dratio3 =0; 

else Dratio3 = DeltaNorm3/Dnorm3; 

if(Dnorm4=0) Dratio4 =0; 

else Dratio4 = DeltaNorm4/Dnorm4; 

if(Dnorm5=0) Dratio5 =0; 

else Dratio5 = DeltaNorm5/Dnorm5; 

if(Unorml > maxUnorm) 

maxUnorm = Unorml; 
else if(Unorm2 > maxUnorm) 
maxUnorm = Unorm2; 
else if(Unorm3 > maxUnorm) 
maxUnorm = Unorm3; 
else if(Unorm4 > maxUnorm) 
maxUnorm = Unorm4; 
else if(Unorm5 > maxUnorm) 
maxUnorm = Unorm5; 


if(Dnorml > maxDnorm) 

maxDnorm = Dnorml; 
else if(Dnorm2 > maxDnorm) 
maxDnorm = Dnorm2; 
else if(Dnorm3 > maxDnorm) 
maxDnorm = Dnorm3; 
else if(Dnorm4 > maxDnorm) 
maxDnorm = Dnorm4; 
else if(Dnorm5 > maxDnorm) 
maxDnorm = Dnorm5; 

inFile_din.getline(buffer,2093); 

inFile_din.getline(buffer,2093); 

inFile_din.getline(buffer,2093); 

inFile_din.getline(buffer,2093); 

inFile_din.getline(buffer,2093); 

inFile_din.close(); 
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// End of all calculations and input file is closed 

// Opening the input file again to add the Noise 

//ifstream inFile_Noise( M NormalizedWhiteNoise.dat ? ’,ios: :in); 
inF ile_din . open(argv [ 1 ] ,ios : : in) ; 


//ofstream outFile_din( M phasel_plus_noise.din M ,ios::out); 
ofstream outFile_din(argv[4],ios: :out); 

inFile_din.getline(buffer,2093); 

outFile_din«setw(5)«buffer«endl; 

inF ile_din»intT mp ; 
outF ile_din«intT mp ; 

inF ile_din»max_node ; 
outFile_din«setw(5)«max_node; 

inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 

inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 

inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 

while(l) 

{ 

inF ile_din»intT mp ; 
outFile_din«setw(5)«intTmp; 

inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 

inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 
if(intTmp == maxnode) 
break; 

} 

while(l) 

{ 

inF ile_din»intT mp ; 
outFile_din«setw(5)«intTmp; 

inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 
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if(intTmp == max_node) 
break; 

} 

inFile_din.getline(buffer,2093); 

outFile_din<<setw(5)<<buffer<<endl; 

inF ile_din»intT mp ; 
outF ile_din«intT mp ; 

inF ile_din»max_element ; 
outFile_din«setw(5)«max_element; 

// int mat; 

inF ile_din»mat ; 

outF ile_din«setw(5)«mat; 

inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 

for(i=0; i<mat; i++) 

{ 

inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 

} 

while(l) 

{ 

inF ile_din»intT mp ; 
outFile_din«setw(8)«intTmp; 
inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 
inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 

if(intTmp = maxelement) 
break; 

} 

inFile_din.getline(buffer,2093); 

outFile_din«setw(5)«buffer«endl; 


// Adding Noise to the UnDamaged Displacements 

ifstream inFile_Noise(argv[2],ios::in); 
if(!inFile_Noise) 

{ 

cerr« M noise file not found! !!!!!! ”«endl; 
return -1; 

} 

done = false; 
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x=0; 


double udNoisel = 0.00; 
double udNoise2 = 0.00; 
double udNoise3 = 0.00; 
double udNoise4 = 0.00; 
double udNoise5 = 0.00; 
double udNoise = 0.00; 

double udSNR, udSNRl, udSNR2, udSNR3, udSNR4, udSNR5; 

while(!done) 

{ 

inF ile_din»intT mp ; 
outFile_din«setw(5)«intTmp; 

if(intTmp == maxnode) 
done = true; 

inF ile_din»intT mp ; 
outFile_din«setw(5)«intTmp; 

inF ile_din»intT mp ; 
outFile_din«setw(5)«intTmp; 

inFile_din»ml ; 
inF ile_N oise»m2 ; 

udNoise = udNoise + m2*m2; 
m= ml *(l+m2*Uratiol *Unorml/maxUnorm); 
outFile_din«setw(25)«setprecision(16)«fixed«m; 

udNoisel = udNoisel + pow(ml*m2*Uratiol*Unorml/maxUnorm,2); 

inF ile_din»intT mp ; 

outFile_din«setw(5)«intTmp; 

inFile_din»ml ; 
inF ile_N oise»m2 ; 

udNoise = udNoise + m2*m2; 

m= ml *(l+m2*Uratio2*Unorm2/maxUnorm); 
outFile_din«setw(25)«setprecision(16)«fixed«m; 

udNoise2 = udNoise2 + pow(ml*m2*Uratio2*Unorm2/maxUnorm,2); 

inF ile_din»intT mp ; 
outFile_din«setw(5)«intTmp; 

inFile_din»ml ; 
inF ile_N oise»m2 ; 

udNoise = udNoise + m2*m2; 
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m= ml *(l+m2*Uratio3*Unorm3/maxUnorm); 
outFile_din<<setw(25)<<setprecision(16)<<fixed<<m; 

udNoise3 = udNoise3 + pow(ml*m2*Uratio3*Unorm3/maxUnorm,2); 

inF ile_din»intT mp ; 
outFile_din«setw(5)«intTmp; 

inFile_din»ml ; 
inF ile_Noise»m2 ; 

udNoise = udNoise + m2*m2; 

m= ml *(l+m2*Uratio4*Unorm4/maxUnorm); 

outF ile_din«setw(25)«setprecision( 1 6)«fixed«m; 

udNoise4 = udNoise4 + pow(ml*m2*Uratio4*Unorm4/maxUnorm,2); 

inF ile_din»intT mp ; 

outFile_din«setw(5)«intTmp; 

inFile_din»ml ; 
inF ile_Noise»m2 ; 

udNoise = udNoise + m2*m2; 

m= ml *(l+m2*Uratio5*Unorm5/maxUnorm); 
outF ile_din«setw(25)«setprecision( 1 6)<<fixed«m«endl; 

udNoise5 = udNoise5 + pow(ml*m2*Uratio5*Unorm5/maxUnorm,2); 

x++; 

} 

udSNRl = 10*logl0(udSignall/udNoisel); 
udSNR2 = 10*logl0(udSignal2/udNoise2); 
udSNR3 = 10*logl0(udSignal3/udNoise3); 
udSNR4 = 10*logl0(udSignal4/udNoise4); 
udSNR5 = 10*logl0(udSignal5/udNoise5); 
udSNR = 10*logl0(udSignal/udNoise); 
inFile_din.getline(buffer,2093); 
outF ile_din«buffer ; 

inFile_din.getline(buffer,2093); 
outF ile_din«setw( 1 )«buffer«endl; 

inFile_din.getline(buffer,2093); 

outFile_din«setw(5)«buffer«endl; 

// Adding Noise to the Damaged Displacements 
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//ifstream inFile_Noise2( ? ’NormalizedPinkNoise.dat ? ’,ios::in); 
ifstream inF ile_N oise2(argv[3 ] ,ios : : in); 

if( ! inF ile_N oise2) 

{ 

cerr« ? ’noise.dat file not found! !!!!!! "«endl; 
return -1; 

} 

done = false; 
x=0; 

double dNoisel = 0.00; 
double dNoise2 = 0.00; 
double dNoise3 = 0.00; 
double dNoise4 = 0.00; 
double dNoise5 = 0.00; 
double dNoise = 0.00; 

double dSNR, dSNRl, dSNR2, dSNR3, dSNR4, dSNR5; 

while(!done) 

{ 

inF ile_din»intT mp ; 
outFile_din«setw(5)«intTmp; 

if(intTmp == maxnode) 
done = true; 

inF ile_din»intT mp ; 
outFile_din«setw(5)«intTmp; 

inF ile_din»intT mp ; 
outFile_din«setw(5)«intTmp; 

inFile_din»ml ; 
inFile_Noise2»m2; 

dNoise = dNoise + m2*m2; 
m= m 1 * ( 1 +m2 * Dratio 1 * Dnorm 1 /maxDnorm) ; 
outFile_din«setw(25)«setprecision(16)«fixed«m; 

dNoisel = dNoisel + pow(ml*m2*Dratiol*Dnorml/maxUnorm,2); 

inF ile_din»intT mp ; 

outFile_din«setw(5)«intTmp; 

inFile_din»ml ; 
inFile_Noise2»m2; 

dNoise = dNoise + m2*m2; 

m= ml *(l+m2*Dratio2*Dnorm2/maxDnorm); 
outFile_din«setw(25)«setprecision(16)«fixed«m; 
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dNoise2 = dNoise2 + pow(ml*m2*Dratio2*Dnorm2/maxUnorm,2); 

inF ile_din»intT mp ; 
outFile_din«setw(5)«intTmp; 

inFile_din»ml ; 
inF ile_Noise2»m2 ; 

dNoise = dNoise + m2*m2; 

m= ml *(l+m2*Dratio3*Dnorm3/maxDnorm); 
outFile_din<<setw(25)<<setprecision(16)<<fixed<<m; 

dNoise3 = dNoise3 + pow(ml*m2*Dratio3*Dnorm3/maxUnorm,2); 

inF ile_din»intT mp ; 

outFile_din«setw(5)«intTmp; 

inFile_din»ml ; 
inF ile_Noise2»m2 ; 

dNoise = dNoise + m2*m2; 

m= ml *(l+m2*Dratio4*Dnorm4/maxDnorm); 
outFile_din«setw(25)«setprecision(16)«fixed«m; 

dNoise4 = dNoise4 + pow(ml*m2*Dratio4*Dnorm4/maxUnorm,2); 

inF ile_din»intT mp ; 
outFile_din«setw(5)«intTmp; 

inFile_din»ml; 
inF ile_Noise2»m2 ; 

dNoise = dNoise + m2*m2; 

m= ml *(l+m2*Dratio5*Dnorm5/maxDnorm); 
outFile_din«setw(25)«setprecision(16)«fixed«m«endl; 

dNoise5 = dNoise5 + pow(ml*m2*Dratio5*Dnorm5/maxUnorm,2); 
x++; 

} 

dSNRl = 10*logl0(dSignall/dNoisel); 
dSNR2 = 10*logl0(dSignal2/dNoise2); 
dSNR3 = 10*logl0(dSignal3/dNoise3); 
dSNR4 = 10*logl0(dSignal4/dNoise4); 
dSNR5 = 10*logl0(dSignal5/dNoise5); 

dSNR= 10*logl0(dSignal/dNoise); 

double SNR, Signal, Noise; 

Signal = udSignal+dSignal; 
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Noise = udNoise+dNoise; 


SNR= 10*logl0(Signal/Noise); 

cout« M \nDISPLACEMENTS TO NOISE RATIO (SNR):”; 
cout«setw( 1 0)«SNR«endl; 

cout«”\nAVERAGE MAGNITUDE (RMS VALUE) OF ” 

«”THE NOISE ADDED TO THE DISPLACEMENTS: 
cout«sqrt(Noise/(l 68 1 *5))«endl; 

inFile_din.getline(buffer,2093); 
outF ile_din«buffer ; 

inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 

inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 

inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 

inFile_din.getline(buffer,2093); 
outF ile_din«buffer«endl; 

inFile_Noise.close(); 
inFile_Noise2.close(); 
inFile_din.close(); 
outF iledin. close() ; 

return 0; 


A.2 Generation of White Noise and Pink Noise 

// GENERATION OF WHITE NOISE AND PINK NOISE 

#include<iostream> 

#include<iomanip> 

#include<cstdlib> 

#include<ctime> 

#include<fstream> 

using namespace std; 

//int main() 
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int main(int argc ,char * argv[]) 

{ 

if(argc !=3 ) 

{ 

cout« ? ’Error 

cout« ? ’Usage: commandline< White_Noise_file Pink_Noise_file >"«endl; 
return -1; 

} 

// Dynamic memory allocation for an Array 
float *white_noisel; 
int SizeOfData; 

cout«"\nENTER THE NUMBER OF NODAL DISPLACEMENTS:”; 
cin»SizeOfData; 

white_noisel = new float[SizeOfData]; 

//End of allocation process 

float *white_noise2; 

white_noise2 = new float[SizeOfData]; 

float *white_noise3; 

white_noise3 = new float[SizeOfData]; 

float *white_noise4; 

white_noise4 = new float[SizeOfData]; 

float *white_noise5; 

white_noise5 = new float[SizeOfData]; 

// Generating the white noise using rand() function 

srand(time (0)); 

for( int i=0; i<SizeOfData; i++) 

{ 

white_noisel[i] = rand(); 
white_noise2[i] = rand(); 
white_noise3[i] = rand(); 
white_noise4[i] = rand(); 
white_noise5[i] = rand(); 

// NORMALIZING THE WHITE NOISE 

int mini, maxi = 0; 
int min2, max2 = 0; 
int min3, max3 = 0; 
int min4, max4 0; 
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int min5, max5 = 0; 


//float normalized_white_noise[SizeOfData]; 
float *normalizedjwhite_noisel; 
normalized_white_noisel = new float[SizeOfData] 

float *normalized_white_noise2; 
normalized_white_noise2 = new float[SizeOfData] 

float *normalized_white_noise3; 
normalized_white_noise3 = new float[SizeOfData] 

float *normalized_white_noise4; 
normalized_white_noise4 = new float[SizeOfData] 

float *normalized_white_noise5; 
normalized_white_noise5 = new float[SizeOfData] 

for(i=0; i<SizeOfData; i++) 

{ 

if(white_noisel[i] > maxi) 

maxi = white_noisel[i]; 

if(white_noise2[i] > max2) 

max2 = white_noise2[i]; 

if(white_noise3[i] > max3) 

max3 = white_noise3[i]; 

if(white_noise4[i] > max4) 

max4 = white_noise4[i]; 


if(white_noise5[i] > max5) 

max5 = white_noise5[i]; 


} 

mini = maxi; 
min2 = max2; 
min3 = max 3; 
min4 = max4; 
min5 = max5; 


for(i=0; i<SizeOfData; i++) 

{ 

if(white_noisel[i] < mini) 

mini = white_noisel[i]; 

if(white_noise2[i] < min2) 

min2 = white_noise2[i]; 

if(white_noise3[i] < min3) 

min3 = white_noise3[i]; 
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} 


if(white_noise4[i] < min4) 

min4 = white_noise4[i]; 

if(white_noise5[i] < min5) 

min5 = white_noise5[i]; 


//Opening the same buffer file to create another file 
//outFile.open( M NormalizedWhiteNoise.daf’,ios::out); 


ofstream outF ile 1 (argv [ 1 ] ,ios : : out) ; 

float factorl, factor2, factor3, factor4, factor5; 
float percent 1 , per cent 2 , percent 3 , percent 4 , percent 5 ; 


float factor; 
float percent; 

int midi, mid2, mid3, mid4, mid5; 

midi = (minl+maxl ) /2 ; 
mid2 = (min2+max2 ) /2 ; 
mid3 = (min3+max3 ) /2 ; 
mid4 = (min4+max4 ) /2 ; 
mid5 = (min5+max5 ) /2 ; 


/* cin>>percent ; 

percent = percent/2; 
f actor=percent/100 ; */ 

cout«" \nENTER THE PERCENTAGE OF THE NOISE LEVEL IN X- 
DIS PLACEMENTS : " ; 

cin>>percentl ; 

cout«"\nENTER THE PERCENTAGE OF THE NOISE LEVEL IN Y- 
DISPLACEMENTS: " ; 

cin>>percent2 ; 

cout«"\nENTER THE PERCENTAGE OF THE NOISE LEVEL IN Z- 
DISPLACEMENTS: "; 

cin>>percent3 ; 

cout«"\nENTER THE PERCENTAGE OF THE NOISE LEVEL IN R1-R0TATI0NS 

n . 
r 

cin>>percent4 ; 

COUt«" \nENTER THE PERCENTAGE OF THE NOISE LEVEL IN R2 -ROTATIONS 

n . 
r 

cin>>percent5 ; 
factorl=percentl/200 ; 
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factor2=percent2/200 ; 
factor3=percent3/200; 
factor4=percent4/200 ; 
factor5=percent5/200 ; 

for(i=0; i<SizeOfData; i++) 

{ 

if(white_noisel [i]<=midl) 

normalized_white_noisel[i] = (white_noisel[i] - mid l)/(midl -mini); 
else 

normalized_white_noisel[i] = (white_noisel[i] - mid l)/(maxl -midi); 


if(white_noise2[i]<=mid2) 

normalized_white_noise2[i] = (white_noise2[i] - mid2)/(mid2-min2); 
else 

normalized_white_noise2[i] = (white_noise2[i] - mid2)/(max2-mid2); 


if(white_noise3 [i]<=mid3) 

normalized_white_noise3[i] = (white_noise3[i] - mid3)/(mid3-min3); 
else 

normalized_white_noise3[i] = (white_noise3 [i] - mid3)/(max3-mid3); 


if(white_noise4[i]<=mid4) 

normalized_white_noise4[i] = (white_noise4[i] - mid4)/(mid4-min4); 
else 

normalized_white_noise4[i] = (white_noise4[i] - mid4)/(max4-mid4); 


if(white_noise5 [i]<=mid5) 

normalized_white_noise5[i] = (white_noise5[i] - mid5)/(mid5-min5); 
else 

normalized_white_noise5[i] = (white_noise5[i] - mid5)/(max5-mid5); 
outFile 1 «setw(20)«setprecision( 1 6)«fixed«normalized_white_noise 1 [i] * factor 1 «" 


«setw(20)«setprecision( 1 6)«fixed«normalized_white_noise2 [i] * factor2«" " 

<<setw(20)<<setprecision(16)<<fixed<<normalized_white_noise3[i]*factor3<<" " 

<<setw(20)<<setprecision(16)<<fixed<<normalized_white_noise4[i]*factor4<<" M 

«setw(20)«setprecision(16)«fixed«normalized_white_noise5[i]*factor5« M " 

«endl; 

} 

outFile l.closeQ; //closing the file 
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// Filtering the White Noise to generate the Pink Noise 

//Assuming that: 

//data is 8 bits 

//sample rate for data is 44 KHz and 
//low pass F is set at 5.5 KHz 
//Therefore 44/5.5 = 8 

#defme F 8 

int avgl, avg2, avg3, avg4, avg5; 
int cnt; 

//float pink_noise[SizeOfData]; 

float *pink_noise 1 ; 

pink_noisel = new float[SizeOfData]; 

float *pink_noise2; 

pink_noise2 = new float [SizeOfData]; 

float *pink_noise3; 

pink_noise3 = new float [SizeOfData]; 

float *pink_noise4; 

pink_noise4 = new float[SizeOfData]; 

float *pink_noise5; 

pink_noise5 = new float[SizeOfData]; 

//Opening the same buffer file to create another file 
//outFile.open("pinknoise.dat",ios::out); 

for (i = 0; i < SizeOfData; i++) 

{ 

for (avgl = 0, avg2 = 0, avg3 = 0, avg4 = 0, avg5 = 0, cnt = 0; 
cnt < F && (i + cnt) < SizeOfData; cnt++) 

{ 

avgl += white_noisel[i + cnt]; //average of several samples 

avg2 += white_noise2[i + cnt]; 

avg3 += white_noise3[i + cnt]; 

avg4 += white_noise4[i + cnt]; 

avg5 += white_noise5[i + cnt]; 

} 

pink_noisel[i] = avgl/cnt; //generation of pink noise 
pink_noise2[i] = avg2/cnt; 

pink_noise3[i] = avg3/cnt; 
pink_noise4[i] = avg4/cnt; 
pink_noise5[i] = avg5/cnt; 


} 
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// NORMALIZING THE PINK NOISE 


//float normalized_pink_noise[SizeOfData]; 
float *normalized_pink_noisel; 
normalized_pink_noisel = new float[SizeOfData]; 

float *normalized_pink_noise2; 
normalized_pink_noise2 = new float[SizeOfData]; 
float *normalized_pink_noise3; 
normalized_pink_noise3 = new float[SizeOfData]; 

float *normalized_pink_noise4; 
normalized_pink_noise4 = new float[SizeOfData]; 

float *normalized_pink_noise5; 
normalized_pink_noise5 = new float [SizeOfData]; 

mini, maxi = 0; 
min2, max2 = 0; 
min3, max3 = 0; 
min4, max4 = 0; 
min5, max5 = 0; 


for(i=0; i<SizeOfData; i++) 

{ 

if(pink_noisel[i] > maxi) 

maxi = pink_noisel[i]; 

if(pink_noise2[i] > max2) 

max2 = pink_noise2[i]; 

if(pink_noise3[i] > max3) 

max3 = pink_noise3[i]; 

if(pink_noise4[i] > max4) 

max4 = pink_noise4[i]; 


if(pink_noise5[i] > max5) 

max5 = pink_noise5[i]; 


} 

mini = maxi; 
min2 = max 2; 
min3 = max 3; 
min4 = max4; 
min5 = max 5; 


for(i=0; i<SizeOfData; i++) 

{ 
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if(pink_noisel[i] < mini) 

mini = pink_noisel[i]; 

if(pink_noise2[i] < min2) 

min2 = pink_noise2[i]; 
if(pink_noise3[i] < min3) 

min3 = pink_noise3[i]; 

if(pink_noise4[i] < min4) 

min4 = pink_noise4[i]; 


if(pink_noise5[i] < min5) 

min5 = pink_noise5[i]; 


} 

midi = (minl+maxl)/2; 
mid2 = (min2+max2)/2; 
mid3 = (min3+max3)/2; 
mid4 = (min4+max4)/2; 
mid5 = (min5+max5)/2; 


//Opening the same buffer file to create another file 
//outFile.open( M NormalizedPinkNoise.daf’,ios::out); 


ofstream outF ile2(argv[2] ,ios : : out) ; 
for(i=0; i<SizeOfData; i++) 

{ 

if(pink_noisel[l]<= midi) 

normalizedjpink_noisel[i] = (pink noiselfi] - mid l)/(midl -mini); 
else 

normalized_pink_noisel[i] = (pink_noisel[i] - mid l)/(maxl -midi); 


if(pink_noise2[l ]<= mid2) 

normalized_pink_noise2[i] = (pink_noise2[i] - mid2)/(mid2-min2); 
else 

normalized jink_noise2[i] = (pink_noise2[i] - mid2)/(max2-mid2); 


if(pink_noise3[l]<= mid3) 

normalized_pink_noise3[i] = (pink_noise3[i] - mid3)/(mid3-min3); 
else 

normalized_pink_noise3[i] = (pink_noise3 [i] - mid3)/(max3-mid3); 


if(pink_noise4[ 1 ]<= mid4) 

normalized_pink_noise4[i] = (pink_noise4[i] - mid4)/(mid4-min4); 
else 

normalized_pink_noise4[i] = (pink_noise4[i] - mid4)/(max4-mid4); 


if(pink_noise5[l]<= mid5) 

normalized jink_noise5[i] = (pink_noise5[i] - mid5)/(mid5-min5); 
else 

normalized_pink_noise5[i] = (pink_noise5[i] - mid5)/(max5-mid5); 
outF ile2«setw(20)«setprecision( 1 6)«fixed«normalized_pink_noise 1 [i] * factor 1 
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«setw(20)«setprecision(16)«fixed«normalized_pink_noise2[i]*factor2« M " 

«setw(20)«setprecision( 1 6)«fixed«normalized_pink_noise3 [i] * factor3 « ,? " 

«setw(20)«setprecision(16)«fixed«normalized_pink_noise4[i]*factor4« M " 

«setw(20)«setprecision(16)«fixed«normalized_pink_noise5[i]*factor5« M " 

«endl; 

} 

outFile2.close(); //closing the file 


delete [] white_noisel; 
delete [] white_noise2; 
delete [] white_noise3; 
delete [] white_noise4; 
delete [] white_noise5; 


delete [] normalized_white_noisel; 
delete [] normalized_white_noise2; 
delete [] normalized_white_noise3; 
delete [] normalized_white_noise4; 
delete [] normalized_white_noise5; 

delete [] pink_noisel; 
delete [] pink_noise2; 
delete [] pink_noise3; 
delete [] pink_noise4; 
delete [] pink_noise5; 


delete [] normalized_pink_noisel; 
delete [] normalized_pink_noise2; 
delete [] normalized_pink_noise3; 
delete [] normalized_pink_noise4; 
delete [] normalized_pink_noise5; 


return 0; 

} 


A.3 Generation of Noised Data Output File and Noised Damage Parameter Signal and also Denoised 
Data Output File Which is Processed using Matlab 

// GENERATION OF NOISED DATA OUTPUT FILE and // 

// DENOISED OUTPUT FILE WHICH IS PROCESSED USING MATLAB // 
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//DESCRIPTION 

//This program adds the noise to the smooth damage parameter values of data output file 
//and stores them in different file to import in Matlab as a signal. 

//And also creates the denoised data output file after denoising the smooth damage //parameter values 
using Matlab Wavelet Toolbox. 

#include<iostream> 

#include<iomanip> 

#include<cstdlib> 

#include<fstream> 

#include<math.h> 

#include<ctime> 

using namespace std; 

//int main() 

int main(int argc ,char * argv[]) 

{ 

int ele, intmp, choice; 

double dtmp, noise, percent, factor, NSDP; 

char buffer[2093]; 

char checkl[] = "DAMAGE PARAMETER FOR M O D E = 1"; 

double maxPai-0.0000; 


while(l) 

{ 


cout«"TO CREATE NOISED DATA OUTPUT FILE : l"«endl«endl 

«"TO CREATE DENOISED DATA OUTPUT FILE : 2"«endl«endl 
«"ENTER YOUR CHOICE : 
cin»choice; 

if(choice !=1 && choice != 2) 

{ 

cout«"\nENTER 1 OR 2 FOR YOUR CHOICE\n"«endl; 
return -1; 

} 

else 

break; 


} 


//CHECKING THE COMMAND LINE ARGUMENTS 
if(choice==l) 


{ 


if(argc != 4) 

{ 

cout«"Error 

cout«"Usage: commandline< dotFile result_dotFile result_NSDPfile>"«endl; 
return -1; 

} 


} 

else 
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{ 

if(argc != 4) 

{ 

cout«"Error 

cout« ? ’Usage: commandline< dotFile Denoised SDPfile 
resultdenoiseddotFile > ?, «endl; 

return -1; 

} 

} 

cout«"\nENTER THE NUMBER OF ELEMENTS = 
cin»ele; 


//ADDING THE NOISE TO THE DAMAGE PARAMETER VALUES AND 
//CREATING NOISED DAMAGED PARAMETER FILE TO DENOISE USING THE //MATLAB 
WAVELET TOOLBOX 

if(choice = 1) 

{ 

//ifstream in_dot(" test, dot", ios::in); 
ifstream in_dot(argv[l], ios::in); 
if(!in_dot) 

{ 

cerr« ?, .dot file not found! !!!!!! M «endl; 

} 


while(l) 

{ 

in_dot.getline(buffer,2093); 

int ch = strcmp(buffer, checkl); 
if(ch==0) 

{ 

in_dot.getline(buffer,2093); 

in_dot.getline(buffer,2093); 

in_dot.getline(buffer,2093); 

in_dot.getline(buffer,2093); 

//RAW DAMAGE PARAMETER 
while(l) 

{ 

in_dot»intmp; 

in_dot»dtmp; 

in_dot»dtmp; 

in_dot»dtmp; 

in_dot»dtmp; 

in_dot»dtmp; 

in_dot»dtmp; 

if(intmp==ele) 
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break; 

} 

indot. getline(buffer,2093 ); 
in_dot.getline(buffer,2093); 
in_dot.getline(buffer,2093); 
indot. getline(buffer,2093 ); 
indot. getline(buffer,2093 ); 


//SMOOTH DAMAGE PARAMETER.... (Finding maximum in the values) 
while(l) 

{ 

in_dot»intmp; 

in_dot»dtmp; 

in_dot»dtmp; 

in_dot»dtmp; 

in_dot»dtmp; 

in_dot»dtmp; 

if(dtmp > maxPar) 

maxPar = dtmp; 


break; 

} 


if(intmp==ele) 

break; 


} 

indot.closeQ; 


in_dot.open(argv[ 1 ],ios: :in); 
if(!in_dot) 

{ 

cerr«" .dot file not found! !!!!!! "«endl; 
return - 1 ; 

} 


srand(time (0)); 

cout«"\nENTER THE PERCENTAGE OF THE NOISE LEVEL = "; 
cin»percent; 

// percent = percent/2; 

factor = percent/ 100; 

// ofstream out_dot("noisetest.dot", ios::out); 

ofstream out_dot(argv[2], ios::out); 

// ofstream out_matlab("NSDPvectors.dat", ios::out); 

ofstream out_matlab(argv[3], ios::out); 

while(l) 

{ 

in_dot.getline(buffer,2093); 
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out_dot«buffer«endl; 

int ch = strcmp(buffer, checkl); 
if(ch==0) 

{ 

in_dot.getline(buffer,2093); 
out_dot«buffer«endl; 
indot. getline(buffer,2093 ); 
out_dot«buffer«endl; 
indot. getline(buffer,2093 ); 
out_dot«buffer«endl; 
in_dot.getline(buffer,2093); 
out_dot«buffer«endl; 

//RAW DAMAGE PARAMETER 
while(l) 

{ 

in_dot»intmp; 

out_dot«setw(6)«intmp; 

in_dot»dtmp; 

out_dot<<setw(20)<<setprecision(12)<<fixed<<dtmp; 

in_dot»dtmp; 

out_dot«setw(20)«setprecision(12)«fixed«dtmp; 

in_dot»dtmp; 

out_dot<<setw(20)<<setprecision(12)<<fixed<<dtmp; 

in_dot»dtmp; 

out_dot«setw(20)«setprecision(12)«fixed«dtmp; 

in_dot»dtmp; 

out_dot<<setw(20)<<setprecision(12)<<fixed<<dtmp; 

//RANDOM NUMBER GENERATION 

noise = 2*(0.5-((double)rand()/(double)(RAND_MAX+l))); 

noise = maxPar*noise*factor; 


} 


in_dot»dtmp; 

out_dot«setw(20)«setprecision( 1 2)«fixed 
«dtmp+noise«endl; 

if(intmp==ele) 

break; 


in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 
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in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

while(l) 

{ 

in_dot»intmp; 
out_dot«setw(6)«intmp ; 

in_dot»dtmp; 

out_dot«setw(20)«setprecision(12)«fixed«dtmp; 

in_dot»dtmp; 

out_dot«setw(20)«setprecision(12)«fixed«dtmp; 

in_dot»dtmp; 

out_dot<<setw(20)<<setprecision(12)<<fixed<<dtmp; 

in_dot»dtmp; 

out_dot<<setw(20)<<setprecision(12)<<fixed<<dtmp; 

//RANDOM NUMBER GENERATION 
noise = 2*(0.5- 

((double)rand()/(double)(RAND_MAX+ 1 ))); 

noise = maxPar*noise*factor; 


} 


break; 


double ndtmp; 
in_dot»dtmp; 
ndtmp = dtmp+noise; 

out_dot«setw(20)«setprecision(9)«fixed 
«ndtmp«endl ; 

out_matlab«setw(20)«setprecision(9)«fixed 

«ndtmp«endl; 

if(intmp=ele) 

break; 


while(l) 

{ 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 
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} 


if(in_dot.eof()) 

break; 


in_dot.close(); 

out_dot.close(); 

out_matlab.close(); 

} 

//CREATING THE DATA OUTPUT FILE WITH DENOISED DAMAGE //PARAMETER 
VALUES 


if(choice == 2) 


// ifstream in_dot("test.dot", ios::in); 

ifstream in_dot(argv[l], ios::in); 
if(!in_dot) 

{ 

cerr«" .dot file not found! !!!!!! "«endl; 

} 

// ifstream in_matlab(”NSDPvectors.dat”, ios::in); 

ifstream in_matlab(argv[2], ios::in); 
if(!in_matlab) 

{ 

cerr« n denoised file not found! !!!!!! "«endl; 

} 

// ofstream out_dot( M Denoisetest.dot", ios::out); 

ofstream out_dot(argv[3], ios::out); 

while(l) 

{ 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

int ch = strcmp(buffer, checkl); 
if(ch==0) 

{ 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

//RAW DAMAGE PARAMETER (denoising only "smooth” for now) 

while(l) 

{ 
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in_dot»intmp; 

out_dot«setw(6)«intmp; 

in_dot»dtmp; 

out_dot<<setw(20)<<setprecision(12)<<fixed<<dtmp; 

in_dot»dtmp; 

out_dot«setw(20)«setprecision(12)«fixed«dtmp; 

in_dot»dtmp; 

out_dot<<setw(20)<<setprecision(12)<<fixed<<dtmp; 

in_dot»dtmp; 

out_dot«setw(20)«setprecision(12)«fixed«dtmp; 

in_dot»dtmp; 

out_dot«setw(20)«setprecision(12)«fixed«dtmp; 

in_dot»dtmp; 

out_dot«setw(20)«setprecision( 1 2)«fixed 
«dtmp«endl; 


if(intmp=ele) 

break; 

} 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

while(l) 

{ 

in_dot»intmp; 
out_dot«setw(6)«intmp ; 

in_dot»dtmp; 

out_dot«setw(20)«setprecision(12)«fixed«dtmp; 

in_dot»dtmp; 

out_dot«setw(20)«setprecision(12)«fixed«dtmp; 

in_dot»dtmp; 

out_dot«setw(20)«setprecision(12)«fixed«dtmp; 

in_dot»dtmp; 

out_dot«setw(20)«setprecision(12)«fixed«dtmp; 
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in_dot»dtmp; 
in_matlab»N SDP; 


out_dot«setw(20)«setprecision(12)«fixed«NSDP«endl; 


} 


} 


break; 


if(intmp=ele) 

break; 


while(l) 

{ 

in_dot.getline(buffer,2093); 

out_dot«buffer«endl; 

if(in_dot.eof()) 

break; 


} 

in_dot.close(); 

out_dot.close(); 

in_matlab.close(); 


} 

return 0; 


A. 4 Generation of Data Input File using the displacements/mode shapes from the ABAQUS Modal 
Analysis 

//////////GENERATION THE DATA INPUT FILE FOR THE DAMAGE CODE/////////// 

/////USING THE DISPLACEMENTS FROM THE ABAQUS MODAL ANALYSIS///// 


//PASS THE UNDAMAGED AND DAMAGED DISPLACEMENTS FILE WITHOUT 
//THE MINIMUM, MAXIMUM AND TOTAL AT THE END OF THE FILE 


#include<iostream> 

#include<iomanip> 

#include<cstdlib> 

#include<fstream> 

#include<math.h> 

using namespace std; 

int getVectors(char s[], char o[], int Ex, int Ey); 
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//int main() 

int main(int argc, char *argv[]) 

{ 

if(argc != 4) 

{ 

cout« ? ’Error "; 

cout« ? ’Usage: commandline< Undamaged_disp_file Damaged_disp_file result_file >"«endl; 
return -1; 

} 


int nodelD = 1 ; 
int maxnode, ch; 
int maxelement; 
double 1, w, t; 
int Xele,Yele; 

ch = nodelD; 

cout«"\nNOTE: Undamaged and Damaged displacements files shouldn't have"«endl 
«" the minimum, maximum and total at the end of the file"«endl«endl; 

ofstream outFile_din(argv[3], ios::out); 

// ofstream outFile_din("temp.din", ios::out); 

cout«"ENTER THE LENGTHEN X-DIRECTION) OF THE PLATE = 

cin»l; 

cout«endl; 

cout«"ENTER THE WIDTH (IN Y-DIRECTION) OF THE PLATE = 

cin»w; 

cout«endl; 

cout«"ENTER THE THICKNESS OF THE PLATE = 

cin»t; 

cout«endl; 

cout«"ENTER THE NUMBER OF ELEMENTS IN X-DIRECTION = 

cin»Xele; 

cout«endl; 

cout«"ENTER THE NUMBER OF ELEMENTS IN Y-DIRECTION = 

cin»Yele; 

cout«endl; 

cout«endl; 

// 1 = 7.00; w = 5.25; t = 0.125; Xele = 40; Yele = 40; 

maxnode = ((Xele+l)*(Yele+l))+(Xele*Yele); 

outFile_din«"Generated using the displacements from ABAQUS: " 

«l«""X"«w« l,,, X"«t« ,t " plate"«" with ”«Xele«"X"«Yele 
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«" mesh"«endl; 

outFile_din« ? ’5 ? ’« M "«max_node«" ?, « ? ’l 1110 0 l"«endl; 
outFile_din«"0 0"<<endl<<"0"<<endl; 

double y,x; 

// double xcor,ycor; 

x = 1/Xele; 
y = w/Yele; 

//GENERATING NODAL COORDINATES FOR THE DATA INPUT FILE 

double yinc, xinc; 
yinc=0.0000000; 

while(l) 

{ 

xinc = 0.0000000; 
while(l) 

{ 

outFile_din«setw( 1 0)<<nodeID«setw(6)«"0"«setw( 15) 

<<setprecision(8)«fixed«xinc«setw( 1 5)«setprecision(8) 
<<fixed«yinc«setw( 1 5)<<setprecision(8)«fixed«- 

<<setw(8)«"0"«setw( 1 5)<<setprecision(8)«fixed«t 
<<setw(8)«"0"«endl; 

outFile_din«setw(3 l)«setprecision(8)«fixed«xinc«setw( 1 5) 

<<setprecision(8)<<fixed«yinc«setw( 1 5)«setprecision(8) 

<<fixed«t/2«endl ; 

nodeID++; 

if((nodeID<max_node) && (nodelD-ch == Xele+1)) 

{ 

eh = nodelD; 
while(l) 

{ 

outFile_din«setw(10)«nodeID«setw(6)«"0"«setw(15) 

«"0.00000000"«setw(15)«"0.00000000"«setw(15)«''0.00000000" 

«setw(8)«"0"«setw(15)<<setprecision(8)«t«setw(8)«"0"«endl; 


outFile_din«setw(31)«"0.00000000"«setw(15)«"0.00000000"«setw(15)«"0.00000000"< 

<endl; 

nodeID++; 

if(nodeID-ch = Xele) 

{ 

eh = nodelD; 
break; 

} 
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} 

if((nodeID % (2*Xele+l))=l || nodelD > max_node) 
break; 
xinc=xinc+x; 

} 

yinc=yinc+y; 
if(nodeID > max node) 
break; 

cout«"THE NODAL COORDINATES ARE GENERATED"«endl«endl; 

outFile_din«setw(8)«" 1 "«" "«"1 1 0 0 0 l"«endl; 

outFile_din«setw(8)«max_node-l<<" "«"1 1 0 0 0 0"«endl; 

nodeID=l; 
ch = nodelD; 

//CRAETING BOUNDARY CONDITIONS FOR THE DATA INPUT FILE 
while(l) 

{ 

outFile_din«setw(8)«nodeID<<" "«" -1 -1 -1 -1 -1 0"«endl; 

if((nodeID-ch) == Xele) 

{ 

nodelD = nodeID+(Xele+l); 
ch = nodelD; 

} 

else 

{ 

nodeID++; 

} 

if(nodeID > max node) 
break; 

} 

cout«"THE BOUNDARY CONDITIONS ARE GENERATED"«endl«endl; 
max_element = Xele*Yele; 

//WRITING THE MATERIAL PROPERTIES TO THE DATA INPUT FILE 

outFile_din«"l 0 0 "«(Xele+l)*(Yele+l)«endl; 

outFile_din«"4"«setw(8)«max_element«" " 

«"2 5321 1000100 0"«endl; 
outFile_din«" 1 7.67E-04"«endl; 

outFile_din«"2.96E+07 0.3 0.0 0.0"«endl; 
outFile_din«"2 7.67E-04"«endl; 

outFile_din«" 1.00936+07 0.3 0.0 0.0"«endl; 

int nodel, node2, node3, node4, midnode; 
nodel=l; 
ch = nodel; 

//GENERATING THE NODAL CONNECTIVITY 
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for(int elelD = 1 ; elelD <= max element; eleID++) 

{ 

outFile_din«setw(8)«eleID<<" 5 1 0"«endl; 

node2 = nodel+1; 
node3 = nodel+2*Xele+l; 
node4 = node3+l; 
midnode = nodel+Xele+1; 

outFile_din<<setw(8)<<nodel<<setw(8)<<node2<<setw(8)<<node4 
<<setw(8)«node3«setw(8)«midnode«" 1 0.0"«endl; 

nodel++; 

if(nodel-ch == Xele) 

{ 

nodel = nodel+(Xele+l); 
ch = nodel; 

} 

} 

cout«"ELEMENTS (NODAL CONNECTIVITY) ARE GENERATED ? ' «endl«endl ; 

ifstream infile2(argv[l], ios::in); 
ifstream infile3(argv[2], ios::in); 

char *out = argv[3]; 
outFile_din.close(); 

char *sl = argv[l]; 

getVectors(sl, out, Xele, Yele); //CALLING THE FUNCTION FOR THE 
// UNDAMAGED DISPLACEMENTS 
//GENERATED BY THE ABAQUS 
outF ile_din. open(argv [3 ] ,ios : : app) ; 
outFile_din«"0"«endl; 
outFile_din.close(); 

char *s2 = argv[2]; 

getVectors(s2, out, Xele, Yele); //CALLING THE FUNCTION FOR THE 
//DAMAGED DISPLACEMENTS 
//GENERATED BY THE ABAQUS 
outF ile_din. open(argv [3 ] ,ios : : app) ; 
outFile_din« M 1 "«endl; 
outFile_din«Yele«setw(8)«Xele«endl; 
outFile_din« ? ’0 ?, «endl; 
outFile_din« M 4 M «setw(8)« , ’4 M «endl; 

infile2.close(); 
infile3.close(); 
outF iledin. close() ; 

return 0; 
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// FUNCTION TO GET DAMAGE FORCE VECTORS 
int getVectors(char s[], char o[], int Ex, int Ey) 

{ 

char buffer[2093]; 
char ch[5]; 
double tmp; 
double freq; 
int i = 0; 
bool eof(); 

char checkl[] = '’Value”; 
char check2[] = "Node”; 

int node = 1 ; 
int check = node; 

ifstream ifp; 
ifp.open(s, ios: :in); 

ofstream ofp; 
ofp.open(o, ios::app); 


while(l) 

{ 


} 

while(l) 

{ 


ifp.getline(buffer, 2093,' '); 
if(strcmp(buffer,checkl) == 0) 

{ 

ifp»ch; 

ifp»freq; 

ofp«setprecision( 1 5)«freq«endl; 
break; 

} 


ifp.getline(buffer, 2093,' '); 
int r = strcmp(buffer, check2); 
if( r = 0) 

{ 

ifp.getline(buffer, 2093); 
ifp.getline(buffer, 2093); 
ifp.getline(buffer, 2093); 

/AVRITING THE DISPLACEMENTS TO THE DATA INPUT FILE 


while(l) 

{ 


ifp»tmp; 

ofp<<setw(8)<<node<<setw(6)<<”5”<<setw(6)<<" 1 ”; 
ifp»tmp; 

ofp«setw(20)«setprecision(15)«fixed«tmp«setw(6)<<2; 
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ifp»tmp; 

ofp«setw(20)«setprecision( 1 5)<<fixed«tmp«setw(6)«3 ; 
ifp»tmp; 

ofp«setw(20)«setprecision( 1 5)<<fixed«tmp«setw(6)«4; 
ifp»tmp; 

ofp«setw(20)«setprecision( 1 5)<<fixed«tmp«setw(6)«5 ; 
ifp»tmp; 


} 

} 

ifp.close(); 
ofp.close(); 
return 0; 

ofp«setw(20)«setprecision( 1 5)<<fixed«tmp«endl; 

if((node-check) = Ex) 

{ 

node = node+(Ex+l); 
check = node; 

} 

else 

{ 

node++; 

} 

//if(ifp.eof()) 

if((node-(Ex+ 1 ))==(((Ex+ 1 ) * (Ey+ 1 ))+(Ex*Ey))) 
break; 
i++; 

} 

break; 


} 
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Appendix B 
Wavelet Transforms 

B.l An example of Decomposition of a Signal in Haar Wavelet 

Signal, /= (4, 6, 10, 12, 8, 6, 5, 5) 

Length of signal, N= 2 3 = 8 

Averages: a 1 =(aj,a 2 , > a Nl2) 


(/1+/2) 

= V2 ’ 
\ 

(/1+/2)] 


, 2 J 

(/ 3 +/ 4 ) 

V 

(/ 3 +/ 4 )^ 


2 J 

(//V-l + /tv) 

= >/2 

4/jv-i + /v) 

>/2 

2 


Fluctuations: d x - (d\ , d 2 , ,d N / 2 ) 

"n-T- J 

, C/ 3 -/ 4 ) fr('(/3-/4)'l 

41 ^ '"l — 2 — J 

iNn ‘ =V2 l 5 J 

Level- 1 Decomposition: / = (a*|<i 1 ) 

/ = (5 V2 , 1 1 V2 , 7 V2 , 5 V2 | -2 / V2 , -2 / V2 , 2 / V2 , 0) 

2 2 1 

Level-2 Decomposition: / = (a \d \ d ) 

/ = ( 16,12 1 -6,2 1 -2 / V2 , -2 / V2 , 2 / V2 , 0) 

Level-N Decomposition: / = (a ,v | d N | | d 2 \d x ) 

Inverse Haar Transform from Level- 1 Decomposition: 
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(a 1 | d l ) = 


( 2 ^ + — d ^ @2 ^2 ®2 — ^2 a N 12 ~^~^ N /2 ^ N /2 — ^ N /2 


42 ’ 72 ’ V2 ’ 42 


4i 


4i 


=f 


(«V) = (4, 6, 10,12, 8, 6, 5, 5) =/ 


Decomposition in the form of Scaling and Wavelet Signals: 
Scaling Signal: 


Vl 


-L-U-o 

k42 42 J 

{ 1 1 

o, o, — = , — = ,0, — , 


V 


'42’ 42' 


\ 

0 

J 


^/2=[o,0, 0,0,-^,-^ 


First Level Averages, a m = (f, V,l 


Wavelet Signal: 


w \ = 


V-U....,o 

42 42 


w \ = fo,o,-^,— ^,0,....,0 

42 42 J 


Kl 2 = 


0,0, 0,0,4 =, — ]= 

42 42 


First Level Fluctuations, d l m = (/, 


Hence, 

Averages, a” = ( /, V” ) and Fluctuations, d^ = ( /, ^ 


\ 

y 


B.2 How to create your own wavelet in wavelet toolbox in Matlab 

The wavelet could be defined by some function. Be sure that the function has integral zero, so that it 
satisfies the requirement of being a wavelet. Equivalently, its Fourier transform should be zero at zero. 
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The wavemngr_ command permits to add new wavelets and wavelet families to the predefined one in the 
GUI. However, before you can use the wavemngr command to add a new wavelet, you must: 

Choose the full name of the wavelet family(fn) 

Choose the short name of the wavelet family(fsn) 

Determine the wavelet type(wt) 

Define the orders of wavelets within the given family(nums) 

Bulit a MAT-file or M-file(file) 

For wavelets without FIR filters: Define the effective support. 

The short name of the wavelet family, fsn, must be a string of four characters or less. For Matlab, 
wavelets can be defined through the definition of the wavelet function as an m-file. To describe a wavelet, 
if you use type 4, write a Matlab code whose first line is function [psi, t] = file(lb, ub, n). The rest of that 
file should be the formula for the wavelet. For example, here is Matlab ’s code for the Mexican-hat 
wavelet (second derivative of a Gaussian function), which has interval of effective support as [-5, 5]. 

Function [ outl, out2] = Mexican(lb, ub, n) 

out2 = linespace(lb, ub, n); %WAVELET SUPPORT 

outl = out2. A 2; 

outl = (2/sqrt(3)*pi A 0.25))*exp(-outl/2).*(l-outl); 

The “WAVEMNGR” command: 

WAVEMNGR Wavelet manager. 

WAVEMNGR is a wavelet manager used to add, delete, restore or read wavelets. 
WAVEMNGR('add',FN, FSN, WT, NUMS, FILE) or 
WAVEMNGR('add',FN, FSN, WT, NUMS, FILE, B) or 
WAVEMNGR('add',FN,FSN,WT, {NUMS,TYPNUMS } ,FILE) or 

WAVEMNGR('add', FN, FSN, WT, {NUMS, TYPNUMS}, FILE, B) adds a new wavelet family. 

FN = family name (string). 

FSN = family short name (string). 

WT defines the wavelet type: 

WT = 1 for orthogonal wavelets. 

WT = 2 for biorthogonal wavelets. 

WT = 3 for wavelet with scale function. 

WT = 4 for wavelet without scale function. 

WT = 5 for complex wavelet without scale function. 

If the wavelet is a single one, NUMS = ". 
examples: mexh, morl. 

If the wavelet is part of a finite family of wavelets, NUMS is a string containing a blank separated 
list of items representing wavelet parameters. 

example: bior, NUMS = T.l 1.3 ... 4.4 5.5 6.8'. 

If the wavelet is part of an infinite family of wavelets, NUMS is a string containing a blank 
separated list of items representing wavelet parameters, terminated by the special sequence **. 
examples: 

db, NUMS = '1 23456789 10 **'. 
shan, NUMS = '1-1.5 1-1 1-0.5 1-0.1 2-3 **' 

In these last two cases, TYPNUMS specifies the wavelet parameter input format: 'integer' or 'real' 
or ’string'; the default value is 'integer'. 

examples: db, TYPNUMS = 'integer' 
bior, TYPNUMS = 'real' 
shan, TYPNUMS = ’string' 
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FILE = MAT-file or M-file name (string). 

B = [lb ub] specifies lower and upper bounds of 
effective support for wavelets of type = 3, 4 or 5. 

WAVEMNGR(’der,N), deletes a wavelet or a wavelet family where N is the wavelet name or the 
family short name. 

WAVEMNGR( ? restore ? ) restores the previous wavelets. asc ASCII-file. 
WAVEMNGR( ? restore ? ,IN2) restores the initial wavelets. asc ASCII-file. 

OUT1 = WAVEMNGR( ? read f ) returns all wavelets family names. 

OUT1 = WAVEMNGR( ? read ? ,IN2) returns all wavelet names. 

OUT1 = WAVEMNGR( ? read_asc ? ) returns all wavelets information retrieved from wavelets. asc 
ASCII-file. 
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